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ABSTRACT 

The intrinsic column density (Ah) distribution of quasars is poorly known. At the high obscuration end of 
the quasar population and for redshifts z < 1, the X-ray spectra can only be reliably characterized using broad¬ 
band measurements which extend to energies above 10 keV. Using the hard X-ray observatory NuSTAR, along 
with archival Chandra and XMM-Newton data, we study the broad-band X-ray spectra of nine optically selected 
(from the SDSS), candidate Compton-thick (Ah > 1-5 x 10 24 cm -2 ) type 2 quasars (CTQS02s); five new 
NuSTAR observations are reported herein, and four have been previously published. The candidate CTQS02s 
lie at z < 0.5, have observed [O III] luminosities in the range 8.4 < log(L[o m ]/£©) < 9.6, and show evidence 
for extreme, Compton-thick absorption when indirect absorption diagnostics are considered. Amongst the nine 
candidate CTQS02s, five are detected by NuSTAR in the high energy (8-24 keV) band: two are weakly detected 
at the ss 3d confidence level and three are strongly detected with sufficient counts for spectral modeling (> 90 
net source counts at 8-24 keV). For these NuSTAR-detected sources direct (i.e.. X-ray spectral) constraints 
on the intrinsic AGN properties are feasible, and we measure column densities « 2.5-1600 times higher and 
intrinsic (unabsorbed) X-ray luminosities ss 10-70 times higher than pre-NuSTAR constraints from Chandra 
and XMM-Newton. Assuming the A7/.S7 A//-detected type 2 quasars are representative of other Compton-thick 
candidates, we make a correction to the Ah distribution for optically selected type 2 quasars as measured 
by Chandra and XMM-Newton for 39 objects. With this approach, we predict a Compton-thick fraction of 
/ct = 36±12 %, although higher fractions (up to 76%) are possible if indirect absorption diagnostics are 
assumed to be reliable. 

Subject headings: galaxies: active - galaxies: nuclei - quasars - X-rays 


1. INTRODUCTION 


Much of the cosmic growth of supermassive black holes is 
thought to occur during a phase of luminous, heavily obscured 
accretion: an obscured quasar ph ase (e.g., Fabian|1999| Gilli 


|et al.||2007l |Treister et al.||2009) . However, our current cen¬ 
sus of obscured quasars appears highly incomplete. While 
unobscured quasars were first discovered over 50 years ago 


(|Schmidt|119631 |Hazard et al.||1963|), it is only in the last 
decade that (radio-quiet) obscured quasars have been discov- 


ered 

in large numbers (e.g., Zakamska et al. 

(2003 Hickox 

et al. 

2007 Reyes et al. 2008 Stern et al. 2D 

12; Assef et al. 

2013 

Donoso et al.||2013). Furthermore, it is only very re- 


absorbing column densities of Ah > 1.5 x 10 24 cm 2 ; here- 
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after CT) quasars have begun to be robustly identified at X-ray 
energies (e.g., Comastri et al.|20fT] Gilli et al.|201 1 Gandhi 


et al.|2014]|Lanzuisi et al.|2015a|l. 

Identifying and characterizing heavily obscured quasars is 
important for various reasons. Firstly, many less luminous 
AGNs in the local Universe appear to be CT (~ 20-30% 
of the total population; e.g., |Risaliti et al.l 119991 [Burlon 
et al.pOlT) . While observational constraints are challenging 


for distant quasars, a significant population of luminous CT 


ground (CXB) spectrum (e.g., Comastri et al.| 1995; 

Gilli et al. 

2007 1 Treister et all2009; Draper & Ballantyne |20 

0 Akylas" 

etaL|2012 Ueda et al. 2014 

). Secondly, while the orientation- 

based unified model (e.g., / 

uitonucci 1993; Urry & Padovani 


obscured and CT AGNs observed in the local Universe, it is 
unclear wh ether a unified model or some ev olutionary sce¬ 
nario (e.g., Sanders et al.|1988[ Hopkins et al.|2008i is more 


appropriate at higher luminosities and redshifts. Indeed, the 
observed dependence of AGN obscuration on luminosity sug- 


gests a departure from the unified model (e.g., Ueda et al. 
2003| |Simpsonl2005||Treister et al.|2010[| lwasawa et al. 2012: 

Asset' etal.|2014f Buchner et al.|20151|Lacy et al.|20l3 1 . The 

above issues can be addressed using X-ray studies which aim 

to measure the column density (A 7 h) distribution and CT frac¬ 
tion of obscured quasars, important components of CXB mod¬ 
els and important tools for understanding AGN models (e.g., 
|Fabian et al.|2009[ Draper & Ballantyne|2010| i. 

X-ray studies of heavily obscured quasars are extremely 
challenging. For instance, to-date very few optically se¬ 
lected obscured quasars (i.e., “type 2” quasars or “QS Q2s”; 
the definition of this term is provided in Section |2.1| i have 
been unambiguously confirmed as CT using broad-band X- 
ray measurements extending to high energies (> 10 keV; 
e.g., Gandhi et al.||2014 >. Including the high-energy data is 
crucial. Firstly, the number of counts is inherently low at 
< 10 keV, due to photoelectric absorption of the X-ray contin¬ 
uum, which restricts the accuracy of X-ray spectral modeling 
and may lead to an underestimate of the absorbing column 
density and intrinsic luminosity. Secondly, important diag¬ 
nostic features can be missed if the observed X-ray energy 
window is narrow. Such features include the photoelectric 
absorption cut-off (e.g., at ss 10 keV for a z = 0.2 AGN ab¬ 
sorbed by Ah = 10 24 cm -2 gas), and features of Compton 
reflection/scattering from cold, dense gas. The latter become 
prominent when CT levels of photoelectric absorption deeply 
suppress the primary continuum, revealing strong Fe Ka flu¬ 
orescent line emission at 6.4 keV and a Compton reflection 
“hump” at > 10 keV (e.g.,[George & Fabian [l~991| l, and may 
arise from an extended structure such as the torus of the uni¬ 


fied model (e.g., Ghisellini et aL|1994 ). 

NuSTAR ( [Harrison et al.|2013) , launched in June 2012, has 
further opened our window on the X-ray spectra of obscured 
AGNs, with sensitivity up to 78.4 keV. As the first orbiting 
observatory to focus high-energy (> 10 keV) X-rays, it pro¬ 
vides a two orders of magnitude improvement in sensitivity 
and over an order of magnitude improvement in angular res¬ 
olution relative to the previous-generation > 10 keV obser¬ 
vatories. Recent studies have demonstrated that, in the case 
of heavily obscured quasars, the most accurate constraints on 
the absorbing column density and intrinsic X-ray luminos¬ 
ity come from a combination of both NuSTAR and XMM- 
Newton/Chandra data, which provide the broadest possible 
energy band pass for X-ray spectral modeling (e.g.,[Luo et al. 


2013l|Balokovic et al,|2014t[Del Moro et al.120141 |Lansbury| 

et al.|20l4HGanclhi et al.|20l4K ~ 

In this paper, we extend the w ork of ] 
hereafter LI4) and |Gandhi et al.[ ( j2014 
NuSTAR to study the high-energy emiss 


Lansbury et al,[ ( 2014[ 

_ hereafter G14), using 

• study the high-energy emission of SDSS-selected 
QS02s which are candidates for being CT (i.e., candidate 
“CTQS02s”). The targets were initially selected based on 


[O III] A5007 line emission (Zakamska et al. 2003 Reyes 


g., Heckman et al. 

2005; |LaMassa et al. 

Hainline et al. 20 

3), and subsequently 


NuSTAR using the low-energy X-ray data available (e.g., |Jia| 
et al. 2013f. L14 looked at an exploratory sample of three z = 


0.41-0.49 candidate CTQS02s: one was weakly detected and 
shown to have a high column density of Ah > 5 x 10 23 cm -2 ; 
the remaining two were undetected but shown to have sup¬ 
pressed X-ray luminosities in the high-energy regime, sugges¬ 
tive of CT absorption. G14 showed the lower redshift object 
SDSS J1034+6001 (also known as Mrk 34; z = 0.05) to have 
a column density and intrinsic power an order of magnitude 
greater than those measured with the pm-NuSTAR X-ray data, 
unambiguously revealing the object to be a CTQS02. 

We present new results for a further five targets, bring¬ 
ing the NuSTAR -observed SDSS-selected candidate CTQS02 
sample to a total size of nine objects. For the brightest 
two sources we model the broad-band X-ray spectra, for 
one weakly detected source we characterize the spectrum us¬ 
ing the X-ray band ratio, and for all targets (including non¬ 
detections) we use the X-ray:mid-IR ratio to infer the intrinsic 
AGN properties. The paper is organised as follows: Section[2] 
details the sample selection; Section[3]describes the X-ray and 
multiwavelength data, along with data reduction and analysis 
procedures; Section [4] presents the results of X-ray spectral 
and multiwavelength analyses; and Section[5]discusses the re¬ 
sults for the full sample of nine Nit ST A /Conserved candidate 
CTQS02s in the context of the parent QS02 population, in¬ 
cluding an estimation of the Ah distribution and CT fraction 
for z < 0.5. The cosmology adopted is (CIm, Ha, h) = (0.27, 
0.73, 0.71). Uncertainties and limits quoted throughout the 
paper correspond to the 90% confidence levels (CL), unless 
otherwise stated. 


2. THE QS02 SAMPLE 
2.1. Definitions 

Quasars are rapidly accreting black holes which emit large 
amounts of radiation, and have luminosities which typically 
place them above the knee of AGN luminosity function. Mul¬ 
tiple thresholds exist in the literature for separating quasars 
from less luminous AGNs (e.g., “Seyferts”). According to the 
classical threshold of [Schmidt & Gre~en| ( j 19831 , quasars are 
those objects with absolute /i-band magnitudes of Mb < 
—23. Thus far we have used the term “obscured” rather 
loosely, since it has different implications depending on the 
wavelength regime in question. In the optical band, objects 
are identified as obscured if they show narrow line emis¬ 
sion without broad (e.g., Ha or H/3) components, a result of 
the central broad line region being hidden from the observer. 
These objects are classed as type 2s, or QS02s if the luminos¬ 
ity is at quasar levels (in type Is the broad line components 
are visible). At X-ray energies, objects are identified as ob¬ 
scured or “absorbed” if their X-ray continua show evidence 
for being absorbed by gas along the line-of-sight, with col- 
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umn densities of Nu > 10 22 cm -2 . The objects in this work 
originate from a sample of optically- identified QS02s ( |Za-| 
kamska et al.|2003] Reyes et al.|2008 1 . Several X-ray studies 


at < 10 keVnave now provided evidence that these optically- 
identified QS02s are also absorbed at X-ray energies, with 
many objects showing indirect evidence for being absorbed 
by column dens ities in excess of -/Vh = 1.5 x 10 24 cm ~ 2 (i.e., 
CT columns; |Vignali et al.||2006 2010 Jia et al ||2013[ ). In 
this paper we look at the direct evidence for CT absorption in 
these optically-identified QS02s, from X-ray analyses which 
incorporate spectral information at > 10 keV. 


2 . 2 . Sample Selection 


When selecting a sample of obscured quasars to observe at 
X-ray energies, it is important to select based on an indicator 
of the intrinsic AGN luminosity such that the sample is unbi¬ 
ased and as representative of the general population as possi¬ 
ble. The [O III] A5007 line, one of the strongest emission lines 
readily visible in the optical, is a suitable choice since such 
emission arises from gas on large (~ 100 pc) scales, mini¬ 
mizing the effect of nu clear obscuration. Reyes et aT] (|2008[ 


hereafter R08; see also jZakamska et al.||2003| ) presented the 
largest sample of [O III] -selected QS02s, consisting of 887 
objects selected from the SDSS. R08 defined quasars as hav¬ 
ing observed (i.e., not corrected for extinction) [O III] lumi¬ 
nosities of L\ o m ] > 2 x 10 8 Lq, and identified the quasars 
as type 2s (i.e., QS02s) following the standard optical defini¬ 
tion. For comparison, the classical absolute magnitude cut of 


matel 

et al. 

y to Lt 0 m i > 3 x 10 °Lq for type 1 sources ( 
2003). Subsequent Chandra and XMM-Ne 

Zakamska 
wton stud- 

ies (e.g., Ptak et al. 2006 

{Vignali et al. 2006, 2010, Jia et al. 

2013 

LaMassa et al. 201 

4) have investigated the soft X-ray 


Lple 


the largest subsample (71 objects) investigated by Jia et al. 
(|2013[ hereafter J13). Figure |Tj shows redshift versusT 
for the R08 and J13 samples. 


qo hi] 


For our study, we select from the J13 sample. In order to 
infer information about the overall optically selected QS02 
population, we desire a parameter space for which the J13 
sample is broadly representative of the R08 sample. As such 
we apply redshift and luminosity cuts of z < 0.5 and L [q m i 
> 2.5 x 10 8 Lq, respectively (see Figure |T|). For these z 
and Lp m ] ranges; (1) the z and L[o m ] distributions of the 
J13 sample and the R08 sample are consistent according to 
the Kolmogorov-Smirnov (KS) test (p = 0.64 and 0.09 for 
z and L[Om], respectively); (2) the majority (74%) of the 
J13 sample are either serendipitous sources in the soft X-ray 
( Chandra and XMM-Newton) data or were targeted based on 
their [O III] properties, and should therefore be relatively un¬ 
biased with respect to the X-ray properties of the R08 sam¬ 
ple. We exclude SDSS J0913+4056 (z = 0.442; L [0 m] 
= 2.1 x 10 1o Lq), since this infrared bright AGN is an extreme 
outlier and has been targeted for NuSTAR separately (D. Far- 
rah et al., in preparation). The above cuts leave 42 QS02s 
from J13, 39 of which are detected at < 10 keV (according to 
J13 and |Vignali et al.|20Q6|p0T0| >. 

From the J13 subsample above, we first targeted an initial 
three candidate CTQS02s at z ~ 0.4-0.5 (this subselection 
is described in L14). Since these three objects were weakly 
or not detected with NuSTAR, for the succeeding targets de¬ 
scribed herein greater consideration was given to the predicted 


to 10 


5 w 9 


10 s 
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Z 

Fig. 1.— Observed (i.e., extinction-uncorrected) [O III] A5007 line lumi¬ 
nosity (L [0 in]) versus redshift ( z ). The R08 sample of SDSS-QS02s is 
shown as grey dots. The J13 sample of Chandra- and XMM-Newton-observed 
objects is indicated in blue, with *+’ and ‘x’ symbols indicating < 10 keV 
detectio ns and non detections, res pectively (according to X-ray analyses in 
J13 and |Vignali et al.12006[|2010) . The dashed lines mark out the parameter 
space used in this work (z < 0.5 and L[q m ] > 2.5 x 10 8 Lq), for which 
the J13 sample is broadly representative of the R08 sample. Our NuSTAR- 
observed subsample of candidate CTQS02s is highlighted by black points, 
with circles marking the five recently observed objects presented in this study, 
diamonds marking the three z ~ 0.4-0.5 objects presented in LI4, and the 
triangle marking the low redshift (z = 0.05) object presented in G14. 

NuSTAR 8-24 keV count ratefMThe predictions were achieved 
by extrapolating from the < 10 keV data, assuming a variety 
of physically motivated torus models which cover a range of 
column densities (10 23 < JVh < 10 25 cm -2 ). To the remain¬ 
der of the J13 subsample above, we applied a cut in observed 
X-ray:[O III] luminosity ratio of ^ 2-10 keV A[0 m] < 1 (a 
conservative threshold for targeting the most obscured can¬ 
didates; see section 4.5 in J13), which leaves 12 CT candi¬ 
dates. From this selection, six objects were observed with 
NuSTAR, with preference being given to the objects with high 
8-24 keV count rate predictions. These include the one object 
presented in G14 and the five presented in this paper, bring¬ 
ing the NuSTAR -observed SDSS-selected candidate CTQS02 
sample to a total size of nine objects. 

In this work we present results for the five recently ob¬ 
served candidate CTQS02s SDSS J0758+3923, 0840+3838, 
1218+4706, 1243-0232 and 1713+5729. For the other four 
previously-studied objects (SDSS J0011+0056, 0056+0032, 
1034+6001 and 1157+6003) the detailed reductions and data 
analyses are presented in L14 and G14. Redshifts and [O III] 
luminosities for the five new objects are listed in Table [T| 
The low-energy (< 10 keV) X-ray spectra have previously 
been characterized by J13, who fit the existing Chandra and 
XMM-Newton data with absorbed power law models. For 
SDSS J1218+4706, the column density constrained by J13 
using this direct (i.e.. X-ray spectral) approach is high, but 
less than CT ( N h = 8 . 0 + 4 ® x 10 23 cm -2 ). In the other four 
cases, the directly constrained column densities are compar¬ 
atively low (7 Vh < 3 x 10 22 cm -2 ). This is in strong dis¬ 
agreement with the extremely low X-ray: [O III] ratios, which 
imply CT absorption. J13 recognised this, and thus used indi- 

1 The 8—24 keV band is the standard hard band defined for the NuSTAR 
extragalactic surveys jAlexander et al.|2013|. 
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TABLE 1 

X-ray Observation Log 


Object Name 
(1) 

z 

(2) 

L { o hi] 

(3) 

NuSTAR Observations 


Soft X-ray Observations 


Observation ID 
(4) 

UT Date 

(5) 

^on 

(6) 

^eff 

(7) 

Observatory 

(8) 

Observation ID 

(9) 

UT Date 

(10) 

t 

(11) 

SDSS J075820.98+392336.0 

0.216 

9.02 

60001131002 

2014:255 

48.3 

41.2 

XMM-Newton 

0305990101 

2006:108 

9.1 









0406740101 

2006:295 

14.2 

SDSS J084041.08+383819.8 

0.313 

8.45 

60001132002 

2014:121 

50.5 

38.4 

XMM-Newton 

0502060201 

2007:289 

19.0 

SDSS J121839.40+470627.7 

0.094 

8.56 

60001135002 

2014:145 

41.8 

34.0 

XMM-Newton 

0203270201 

2004:153 

40.8 









0400560301 

2006:321 

43.2 

SDSS J124337.34—023200.2 

0.281 

8.88 

60001136002 

2014:211 

55.5 

46.0 

Chandra 

6805 

2006:115 

10.0 

SDSS J171350.32+572954.9 

0.113 

8.95 

60001137002 

2014:120 

54.5 

45.3 

XMM-Newton 

0305750401 

2005:174 

4.4 


NOTE. - (1): Full SDSS object name. (2): Redshift. (3): Gaussian fit [O III] A5007 line luminosity [log(L[ 0 m]/-k©)L as reported in R08. (4) and (5): NuSTAR 
observation ID and start date (YYYYDDD), respectively. (6): Total on-source time (ks). (7): Effective on-axis exposure time (ks). This is the net value for the 
3-24 keV band, and at the celestial coordinates of the target, after data cleaning. We have accounted for vignetting; despite the sources being “on-axis”, there 
is a small loss of exposure due to the natural dither of the observatory. (8), (9) and (10): Soft X-ray observatory with available data, corresponding observation 
ID(s) and start date(s) (YYYYDDD), respectively. (11): Net on-axis, flaring-corrected exposure time(s) (ks). For XMM-Newton, the quoted value corresponds 
to the EPIC detector used with the longest net exposure time. 


rect diagnostics to estimate the absorption levels. The low TVh 
measurements from direct spectral fitting can be explained as 
due to a combination of the limited energy ranges of Chandra 
and XMM-Newton, low source counts, a nd (es pecially in the 
case of SDSS J1713+5729; see Section [4. 1.3| for further de¬ 
tails) strong contamination at lower energies from other pro¬ 
cesses such as star formation, AGN photoionization, or scat¬ 
tered AGN emission. In the Appendix we give individual ob¬ 
ject information for the five candidate CTQS02s presented in 
this paper, including relevant multiwavelength properties and 
indicators of heavy absorption. In addition, we comment on 
the single NuSTAR-detected candidate CTQS02 from the ex¬ 
ploratory study of L14 (SDSS J0011+0056), for which a close 
inspection of the soft X-ray data reveals strong Fe Ka emis¬ 
sion. 


3. DATA 


This section details the pointed NuSTAR observations and 
data analysis procedures for the five n ewly observed SDSS- 
selected candidate CTQS02s (Section which bring the 
NuSTA /(-observed sample to a total of nine such objects. We 
also deta il the archival Chandra and XMM-Newton data (Sec¬ 
tion |3 .2| ), which facilitate a broad-band X-ray analysis when 
combined with the NuSTAR data. In addition, near-UV to 
mid-IR data from large-area surveys are presented in order 
to characterize the spectral energy distributions (SEDs) of the 
objects and dise ntang le AGN and host galaxy emission in the 
mid-IR (Section [T3| . 


3.1. NuSTAR Data 


The NuSTAR observatory is sensitive at 3-78.4 keV ( |Harri- 
son et a l.|2013) l. The combination of the instrumental back- 
ground and decrease in effective area with increasing energy 
means that 3-«24 keV is the most useful energy band for faint 
sources. NuSTAR consists of two telescopes (A and B), identi¬ 
cal in design, the respective focal plane modules of which are 
referred to as FPMA and FPMB. The point-spread function 
(PSF) has a tight “core” of FWHM = 18" and a half-power 
diameter of 58". 

Table [T] provides details, including dates and exposure 
times, for the most recent five NuSTAR observations of SDSS- 
selected candidate CTQS02s. The data were processed as for 
the L14 sample, using the NuSTAR Data Analysis Software 
(NuSTARDAS) version 1.3.0. For the detected sources, the 


NUPRODUCTS task was used to extract spectra and response 
files. Followi ng other recent N uSTAR studies (Alexa nder et ak] 
|2013| L14; |Luo et al.||2014), we perform photometry in the 
3-24 keV, 3-8 keV, and 8-24 keV bands. The photometry 
is performed for each FPM separately and also for combined 
FPMA+FPMB data (referred to hereafter as “FPMA+B”), to 
increase sensitivity. For source detection, we use prior knowl¬ 
edge of the SDSS coordinates and calculate no-source proba¬ 
bilities assuming binomial statistics (Pb), defining non detec¬ 
tions as Pb > 1% (i.e., < 2.6er). For non detections we calcu¬ 
late upper limits on the net source counts using the Bayesian 
approach outlined in |Kraft et al.| 1991!). For a detailed de¬ 
scription of the source detection and aperture photometry pro¬ 
cedures, we refer the reader to L14. 

Table [2] summarizes the NuSTAR photometry. Two of 
the quasars, SDSS J1218+4706 and 1243-0232, are strongly 
detected; the net source counts for FPMA+B in the 8- 
24 keV band are 188 and 90, respectively. Figure [2] shows 
the 8-24 keV no-source probabilities for the three fainter 
sources, SDSS J0758+3923, 0840+3838 and 1713+5729. 
Poisson, rather than binomial, no-source probabilities have 
been adopted for the purposes of the figure only, to aid inter¬ 
object comparison; these provide a good approximation of 
the binomial no-source pr obab ilitie s (Pb ) since the back¬ 
ground counts are large (Weisskopf et al. 2007). Although 
SDSS J0758+3923 is formally undetected at 8-24 keV, it is 
only just below the adopted detection threshold for this band 
and is weakly detected in the broader 3-24 keV energy band, 
but for FPMA only (P B = 0.63%). SDSS J0840+3838 is 
a non detection. SDSS J1713+5729 is weakly detected with 
FPMA+B for the 8-24 keV band only (P B = 0.22%). In 
general, the detected sources have more net source counts in 
the 8-24 keV band, where the focusing soft X-ray observato¬ 
ries (e.g., Chandra and XMM-Newton ) have little to no sen¬ 
sitivity, than in the 3-8 keV band, where NuSTAR and the 
soft X-ray observatories overlap. This can occur for heav¬ 
ily obscured AGNs, which have extremely flat X-ray spec¬ 
tra and are therefore brighter at > 8 keV. Indeed, the sin¬ 
gle candidate CTQS02 to be detected with NuSTAR in L14, 
SDSS J0011+0056, was only detected in the 8-24 keV band. 
NuSTAR FPMA+B 8-24 keV image cutouts for the three new 
targets detected in this energy band are shown in Figure [3] 
None of these three sources are detected in the most sensitive 
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TABLE 2 

X-ray Photometry: NuSTAR Counts 


Object Name 

Net Counts (3-24 keV) 

Net Counts (3-8 keV) 

Net Counts (8-24 keV) 

SDSS J 

FPMA 

FPMB 

FPMA+B 

FPMA 

FPMB 

FPMA+B 

FPMA 

FPMB 

FPMA+B 

0758+3923 

0840+3838 

1218+4706 

1243-0232 

1713+5729 

SO 4+ 17 - 9 

OU.+-164 

< 25.2 
122.9±f 9 ; 8 3 
56.811®;® 

< 43.1 

< 14.8 

< 17.1 

1 97 0+21.6 
iZ '- Z _20.2 

60 4+ 21 - 7 
OU.4:_20.2 

< 33.5 

< 43.8 

< 28.4 
249.91®®;® 
116. 9 I 2 8 ;® 

< 67.4 

< 29.3 

< 14.6 

32 4+ 12 - 6 

< 32.4 

< 18.1 

< 7.2 

< 8.8 

39 4+!3.3 

-11.8 

< 31.8 

< 13.3 

< 18.1 

< 13.4 

64 7 + 17 - 8 
—16.4 
00 0 +I 8.8 

oo.o _ 17 3 

< 21.5 

< 30.4 

< 19.1 

91 4 +171 

15.6 

40 0 ~*~ 15 ' 8 

^ u - u _ 14 .3 

< 33.9 

< 21.8 

< 21.5 

96 7 + 17 - 7 
aD - ‘ -16.2 

49 6^ 17 ' 2 

+9.0-15 7 

< 36.3 

< 45.0 

< 31.5 

188.0l?li 

89.6 +2 i'3 

38.iia:! 


NOTE. - NuSTAR net source counts for the candidate CTQS02s. FPMA and FPMB are the individual focal plane modules belonging to the two telescopes which 
comprise NuSTAR. “FPMA+B” refers to the combined FPMA+FPMB data. 



Fig. 2.— NuSTAR photometry at 8-24 keV for the faintest three sources, 
SDSS J0758+3923, 0840+3838 and 1713+5729 (circles, squares and dia¬ 
monds, respectively). Gross source counts and background counts (scaled 
to the source aperture) are shown. The dashed lines indicate tracks of con¬ 
stant Poisson no-source probability (a good approximation of P \>. given 
the large background counts considered here; [Weisskopf et al J2007| . The 
solid black line shows our adopted detection threshold of Fg = 1%- 
Only SDSS J1713+5729 is detected: while it is not detected in the in¬ 
dividual FPMs, the increased sensitivity in FPMA+B (i.e., the combined 
FPMA+FPMB data) results in a significant detection, with Pg = 0.22%. 


Swift BAT all-sky catalogs (e.g., Baumgartner et al.||2013j>, 
and direct examination of the 104 month Swift BAT maps 
shows no excess above 2a (for details of the maps and pro¬ 
cedures, see Koss et al. |2013) >. Therefore, NuSTAR has pro¬ 
vided the first real detections of these targets at high energies 
(> 10 keV). 

For the AuS'TAP-detected sources, it is important to rule 
out confusion with and contamination from other nearby 
X-ray sources. Both of these are extremely unlikely: in 
the soft X-ray (Chandra and XMM-Newton ) imaging of the 
M/57i4P-detected sources, the only neighbouring source de¬ 
tected within 88" (i.e., the radial distance containing an 
encircled-energy fraction of ~ 85% for the NuSTAR PSF) of 
the SDSS positions lies at an angular separation of 51" from 
SDSS J1218+4706 (i.e., outside our adopted source aperture 
radius) and is a factor of sa 20 fainter in the XMM-Newton 
energy band. 

Table [3] lists the aperture-corrected NuSTAR fluxes and rest- 
frame 10-40 keV luminosities (L 10-40 kevi uncorrected for 
absorption). The fluxes were obtained using photometry, as¬ 
suming an effective photon index (i.e., for an unabsorbed 
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Fig. 3.— NuSTAR 8-24 keV images for the three objects detected in this 
energy band: SDSS J1218+4706, 1243-0232 and 1713+5729 (left to right, 
respectively). Top row: unsmoothed photon images. Bottom row: images 
smoothed with a top hat function of radius 14 pixels, corresponding to 34" 5 
(for aesthetic purposes only). The 45" radius source apertures are shown, 
centered on the SDSS positions. The major tickmarks indicate 1 arcmin off¬ 
sets in right ascension (R.A.; horizontal axis) and declination (Deck; vertical 
axis). 


power law model) of I ’ e ff = 0.3 and using count rate to 
flux conversion factors which account for the NuSTAR re¬ 
sponse and effective area. Often r e g- = 1.8 (a typical value 


for the 3-24 keV emission of AGNs; e.g., Alexander et al. 
2013) is assumed for such extrapolations, but the NuSTAR- 
detected candidate CTQS02s have extremely flat observed 
spectral slopes at 3-24 keV (see Section[4]i, in agreement with 
r eff = 0.3 in all cases. For each object our measured NuSTAR 
flux is in agreement with the soft X-ray observatory (Chan¬ 
dra or XMM-Newton ) measurement at 3-8 keV, the energy 
band for which the observatories overlap. For the three faint 
or undetected sources (SDSS J0758+3923, 0840+3838 and 
1713+5729), the L 10-40 keV values were obtained by extrap¬ 
olating from the observed-frame 8-24 keV fluxes assuming 
r e ff = 0.3. For the two sources with good NuSTAR photon 
statistics (SDSS J1218+4706 and 1243-0232) the L 10 _ 4 o keV 
values were calculated using the best-fitting spectral models 
(Section [44) . 


3.2. Lower Energy X-ray Data 

To incorporate lower energy (< 10 keV; or “soft”) X-ray 
data in our study, we use archival Chandra and XMM-Newton 
observations, limiting the analysis to the 0.5-8 keV and 0.5- 
10 keV bands, respectively. Table IT] provides details of the 
archival soft X-ray observations, including dates and net ex¬ 
posure times. For the sources with poor photon statistics, we 
perform photometry using identical p roce dures to those for 
the NuSTAR photometry (see Section ET For the sources 
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TABLE 3 

Multiwavelength Flux and Luminosity Measurements 


Object 

Observed-frame Flux (10" 

_1 Q _ 1 

° erg s 

cm -2 ) 

Rest-frame Luminosity (10 42 

erg s 

a 



Chandra / XMM 


NuSTAR 


Chandra / XMM 

NuSTAR 


SED Modeling 


SDSS J 

3-8 keV 

3-8 keV 

3-24 keV 

8-24 keV 

2-10 keV 

10-40 keV 

6 fim 

0.1-30 fim 

6 fim 

( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

(8) 

(9) 

( 10 ) 

0758+3923 

0 13+ 003 
u -- lo - 0.02 

< 0.12 

< 0.69 

< 0.93 

O qq+0.40 
z -°°—0.35 

< 23.22 

347 ± 19 

0.88 ± 0.01 

0.97 ±0.00 

0840+3838 

1218+4706 

1243-0232 

1713+5729 

< 0.13 

U - 0 '-0.47 

0 15+ 008 
u ' lo -0.09 

< 0.30 

< 0.09 

0 49+ 013 
u, ^ y - 0.12 

0 19+ 011 
u - iy - 0.10 

< 0.12 

< 0.48 

4 66" 1-0 ’ 55 

^ oo -0.52 

± ' oo -0.39 

< 0.95 

< 0.69 

4 49 + 0,58 
^•^-0.54 

1 62 +0 ' 41 
1-0 —0.39 

0 69 +0 ' 35 
u ‘° -0.33 

< 3.93 
i-38±S;lg 

5-74 t° 0 f 6 

< 1.07 

< 35.23 
14.00l 44 1 ® 3 
54-60 tlf 7 

*™±Im 

130 ± 10 
73 ±3 

25 ±4 

305 ± 21 

0.63 ± 0.03 
0.91 ± 0.01 

0.30 ± 0.04 

0.92 ± 0.02 

0.73 ± 0.04 
0.99 ±0.01 

0.52 ± 0.07 

n qq+ 0.01 
u ' aM —0.03 


NOTE. - Columns (2) to (7): Hard X-ray ( NuSTAR) and soft X-ray (Chandra or XMM-Newton ) fluxes and luminosities. The rest-frame X-ray luminosities are 
observed values, i.e. uncorrected for absorption, and are in units of 10 42 erg 1 . The NuSTAR fluxes are from photometry in t hree observed-frame energy 
bands, assuming T e fj = 0.3. The rest-frame 10-40 keV luminosities are determined from the best-fitting spectral models (Scction |4. 1 [ for SDSS J1218+4706 
and 1243-0232, and by extrapolating from the observed-frame 8-24 keV band (assuming r e ff = 0.3) for SDSS J0758+3923, 0840+3838 and 1713+5729. The 
Chandra and XMM-Newton fluxes and luminosities are determined from spectroscopy for SDSS J0758+3923, 1218+4706 and 1243-0232, and from aperture 
photometry in the observed-frame 3-8 keV and rest-frame 2-10 keV bands for SDSS J0840+3838 and 1713+5729 (assuming r e fj = 0.3). Columns (8) to 
(10): Best-fit parameters from the near-UV to mid-IR SED modeling in Section [XT] The errors shown correspond to standard deviations from a Monte Carlo 
re-sampling of the photometric data. Column (8): rest-frame 6 /im luminosity for the AGN only, L6/rm (uL v ), in units of 10 42 erg s 1 . This value is intrinsic 
(i.e., corrected for dust extinction). Column (9): The fractional contribution of the AGN to the total integrated intrinsic luminosity between 0.1 and 30 ft, m. 
Column (10): The fractional contribution of the AGN to the observed (i.e., uncorrected for dust extinction) monochromatic rest-frame 6 //m flux. 


with good photon s tatist ics, we model the X-ray s pect ra with 
XsPEC (see Section |4T| >. As mentioned in Section [3TT| source 
confusion is extremely unlikely: there are no neighbouring 
sources detected within 51" of the QS02 positions. Measure¬ 
ments of the observed-frame 3-8 keV fluxes and rest-frame 
2-10 keV luminosities (uncorrected for absorption) are listed 
in Tabled 

For the source with Chandra coverage (SDSS J1243- 
0232), we process the data using CHANDRA.REPRoJ^] The 
source events are extracted from a circular 2 / /5 radius aper¬ 
ture. The background events are extracted from a background 
source-free annulus centered on the source coordinates, with 
an inner radius of 8" and an outer radius of 80". Since 
SDSS J1243-0232 is on-axis, a large fraction (> 90%) of the 
source counts lie within the source aperture. Given this, and 
the extremely low net source counts measured (9), contamina¬ 
tion of the background region by source counts is negligible. 

For the sources with XMM-Newton coverage, we analyze 
data products from the Pipeline Processing System (PPS) us¬ 
ing the Science Analysis Software (SAS v.13.5.0). To deter¬ 
mine appropriate count rate thresholds for background flare 
subtraction, we visually examine the light curves. In all cases 
the fraction of exposure time removed is < 30%, except in the 
case of obsID 0305750401 where the fraction is 49%. The 
exposure times in Table [T| are flaring-corrected. The source 
events are extracted from circular regions of 8-20" radius (de¬ 
pending on source brightness and off-axis angle). The back¬ 
ground events are extracted from regions of area 70 x 70" to 
140 x 140", using either an annulus centered on the source po¬ 
sition or an offset region if it is necessary to avoid chip-gaps 
or nearby sources. We combine the MOS1 and MOS2 data 
using the SAS task EPICSPECCOMBINE, and simultaneously 
fit the PN and MOS data when performing spectral analyses. 

In the case of SDSS J1218+4706, we use the two 
archival XMM-Newton observations with the longest expo¬ 
sures and most recent start dates (obsIDs 0203270201 and 
0400560301). For obsID 0203270201, SDSS J1218+4706 
lies close to the on-axis position. In this instance we only 

2 http://cxc.harvard.edu/ciao/ahelp/chandra_repro.html 


use the MOS data, since the source lies on a chip-gap for PN. 
For obsid 0400560301, SDSS J1218+4706 lies far off-axis. 
In this case we only use the PN data, since the source lies on 
a chip-edge in MOS1, and there are relatively low net counts 
with MOS2 (65). 

3.3. Near-UV to Mid-IR SED Analysis 

Here we analyse near-UV to mid-IR (0.3-30 /im) spectral 
energy distributions (SEDs) for the five candidate CTQS02s 
presented in this work, and the one presented in G14 
(SDSS J1034+6001), with the primary aim of reliably mea¬ 
suring the AGN emission at mid-IR wavelengths. The pho¬ 
tometric data (shown in Figure [4| are collated from the 
SDSS (Data Release 7; |Yorket al |2000j >, the WISE All-Sky 
source catal og ( | Wright et al.|2010] >, and the Spitzer (Werner] 
|et al.|[2004|) Enhanced Imaging Products Source List (for 
SDSS JI243-0232 only). The SDSS fluxes are corrected for 
Galactic extinction. The photometric data adopted are pro¬ 
vided in the Appendix. In order to provide a consistent SED 
analysis across the full sample of nine NuSTAR -observed can¬ 
didate CTQS02s, we use the same SED decomposition pro 
cedure as that applied in L14 to the initial three objects. Fol 


lowing the methodology detailed in Assef et al. (2008 2010 
2013), each SED is modeled as the best-fit, n on-neg ative, 


linear combination of four empirical templates (Assef et al. 
2010| , including one AGN template and three galaxy tem¬ 


plates for: an old stellar population (“elliptical” or E), ongo¬ 
ing star formation (“spiral” or Sbc), and a starburst popula¬ 
tion (“irregular” or Im). The internal dust extinction of the 
AGN component is included as a free parameter in the mod¬ 
eling. The model solutions are shown in Figure 0 and the 
following best-fitting parameters are listed in Table 3] a, the 
fractional contribution of the AGN to the total intrinsic (i.e., 
corrected for the dust extinction of the AGN component) in¬ 
tegrated 0.1-30 pm luminosity; the fractional contribu¬ 
tion of the AGN to the total observed (i.e., uncorrected for the 
dust extinction of the AGN component) monochromatic rest- 
frame 6 pm flux; and the intrinsic AGN luminosity at 

rest-frame 6 pm (vL v ). The errors represent standard devia¬ 
tions from a Monte Carlo re-sampling of the photometric data 
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Fig. 4.— Near-UV to mid-IR spectral energy distributions (SEDs) for the five candidate CTQS02s presented in this work, and the one (SDSS J1034+6001) 
presented in G14. AGN (blue dashed curve) and galaxy (dotted curves) templates were combined in the best-fit modeling of the photometric data (colored circles 
for the sources detected at 8-24 keV with NuSTAR, and grey circles for the 8-24 keV non-detections), following|Assef et al.|(2008||2010||2013ji. The three galaxy 
templates correspond to an old stellar population (“elliptical” or E; red), ongoing star formation (“spiral” or Sbc; purple), and a starburst population (“irregular” 
or Im; orange). The gray curve shows the combined model solution. The systems are all AGN-dominated in the mid-IR waveband based on this analysis, except 
for SDSS J1243-0232, which has comparable contributions from the AGN and the host galaxy; see Table[3] 


over 1000 iterations, and thus account for possible model de¬ 
generacies. In all cases the integrated light properties (i.e., 
the total galaxy and AGN contributions) are well constrained, 
which is required to accurately determine a, fe^m and Lq ^ m . 
Since the primary goal of the SED modeling was to reliably 
measure these parameters, we do not make inferences about 
the host galaxy properties from the best-fit combination of 
hostgalaxy templates. SDSS J1034+6001, not shown in Ta¬ 
ble [3]since the X-ray analysis is presented in G14, has / J(i;iln 

= (1.20 ± 0.09) x 10 44 erg s" 1 , a = 0.90 ± 0.02, and / 6/jrn 

— n qc+o.02 

The a constraints demonstrate that the candidate CTQS02s 
in Figure [4] require an AGN component at a very high confi¬ 
dence level, and that in general the AGN contributes strongly 
to the intrinsic emission across the broad 0.1-30 pm wave¬ 
length range (all but one object have a > 0.6). The high 
/ 6(lm values (all but one have / 6/im > 0.7) indicate that the 
observed monochromatic 6 pm fluxes are AGN-dominated. 
The presence of an AGN at mid-IR wavelengths may also 
be inferred using WISE color diagnostics. In Figure [5] we 
show the six objects from Figure [4} and the three from F14, 
on the WISE W1-W2 (i.e., [3.4 pm]-[4.6 pm]) versus W2- 
W3 (i.e., [4.6 pm]-[12.0 pm]) plane. Generally, sources with 
larger W1—W2 values have stronger AGN contributions. We 
compare with the AGN ‘wedge’ of Mateos et al.] (j2013)> and 
the W1-W2 color cut of Stern et al. (2012|, which may be 


used to identify AGN-dominated systems. Out of the total 
sample of nine candidate CTQS02s, five are AGN-dominated 
according t o both criteria, and one (SDSS J0056+0032 ) falls 
s et al7|(|2013Jl wedge but lies above the Stern 
. This is in good agreement with the SED 


below the 


et al. (20 


Mateos i 
2) cut. 


good agreement 
modeling for these sources, where a > 0.9 in all cases. The 
remaining three sources (SDSS J0011+0056, 0840+3838 and 
1243-0232) fall below both of the selection regions, although 
SDSS J0840+3838 is consistent with satisfying the Stern et ah 
(|2012j) AGN selection criterion given the errors. This sup¬ 
ports the SED modeling, from which it is concluded that these 
three sources are the least AGN dominated (a ss 0.3-0.6, 
and f6fj.ni ~ 0.5-0.7). The WISE colors of the objects agree 
with the expectations; in general, the CTQS02 population ap¬ 
pears to follow the WISE color distribution of the total QS02 
population, with a fraction of objects (~ 70%) lying within 
the AGN wedge ( [Mateos et al.|2013 i. In the local Universe, 
~ 40% of the currently known bona fide CT AGNs lie within 
the wedge (Gandhi et al.|2015 


In addition to the near-UV to mid-IR SED, one of the can¬ 
didate CTQS02s presented in this work (SDSS J1713+5729) 
has a detection at far-IR wavelengths with IRAS which allows 
us to assess the extent to which star formation could con¬ 


tribute to the soft X-ray emission (Section 4.1.3 i. 


4. RESULTS 
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Fig. 5.— WISE color-color diagram for the NuSTAR-observed SDSS- 
selected candidate CTQS02s from this study (circles), G14 (triangle) and 
L14 (diamonds). We compare with the AGN c olor cut o f Stern et al. < |2012[ 
W1—W2 > 0.8) and the AGN ‘wedge’ of |Mateos et al.| (TfFH . The 
filled and empty symbols mark sources which are strongly AGNTIominated 
(d > 0.9) and less AGN-dominated (d < 0.6), respectively, at mid-IR wave¬ 
lengths according to our SED modeling. For the five objects that lie within 
the AGN wedge, the error bars are smaller than the symbols. 



Fig. 6.— Observed X-ray properties of the NuSTAR-obseived candidate 
CTQS02 sample. Left panel: properties measured at low energies with 
Chandra and XMM-Newton. Right panel: properties measured at high ener¬ 
gies with NnSTAR. Only detected sources are shown. Circles, diamonds and 
the triangle indicate the objects presented in this work, L14, and G14, respec¬ 
tively. The effective photon index (r e ff), which provides a basic description 
of the overall X-ray spectral shape, was obtained by fitting an unabsorbed 
power law model to the data for each source. The rest-frame X-ray luminosi¬ 
ties (L^ bs ; bottom axis) are observed values, i.e. uncorrected for absorption. 
For the five objects presented in this paper, the luminosities correspond to 
those in Table p] For the L14 and G14 objects the luminosities have been 
calculated using the same methodology: spectral modelin g wh ere possible, 
or photometry following the procedure outlined in Scction |TT| In the cases 
of the X-ray faint sources SDSS J0011+0056 and SDSS J1713+5729, the 
r efT fo r 3-2 4 keV was estimated from the NuSTAR band ratio (ERxj,,; see 
Section |4~2) . 


To summarise the NuSTAR source detection for the five 
SDSS-selected candidate CTQS02s presented in this work: 
two are strongly detected, one is weakly detected, and two are 
undetected by N uSTAR in the high energy band (8-24 keV). 
In Section 14. 1 we present the results of X-ray spectral fit¬ 
ting with ASPEC for the three brightest objects. In Section 
4.2 we present the X-ray band ratios of all of the NuSTAR- 
detected candidate CTQS02s, comparing to model predic¬ 
tions. For the weakly detected source SDSS J1713+5729, this 
is an appropriate method for characterizing the broad-band 
X-ray spectrum. These two sections give direct (i.e.. X-ray 
spectral ) co nstraints on absorbing column densities (TVh). In 
Section 4.3 we present indirect constraints from a multiwave¬ 
length diagnostic for the entire sample, including NuSTAR 
non-detections. 

First we take a brief look at the overall X-ray spectral 
shapes for the full sample of nine NuSTA /^-observed candi¬ 
date CTQS02s. Figure [6] shows the effective photon indices 
(r eff ), measured through unabsorbed power law fits to the in¬ 
dividual Chandra or XMM-Newton (0.5-8 keV) and NuSTAR 
(3-24 keV) spectra. The spectral shapes observed by Chandra 
and XMM-Newton vary significantly over an order of magni¬ 
tude in (non absorption corrected) rest-frame 2-10 keV lumi¬ 
nosity. The increase in I’ e ff (0.5-8 keV) toward lower lumi¬ 
nosities may reflect an increase in the relative contribution to 
the low-energy spectra from processes unrelated to the direct 
AGN emission, such as thermal plasma emission due to star 
formation or AGN-powered photoionization. In contrast, the 
spectra seen by NuSTAR are consistent with having the same 
effective photon index: excluding upper limits, the mean is 
r eff (3-24 keV) « 0.3. 

4.1. X-ray Spectral Analysis: Best-jit Modeling 

Here we use broad-band X-ray spectral modeling for the 
two brightest NuSTAR-detected sources presented in this pa¬ 
per (SDSS J1218+4706 and SDSS J1243-0232) to measure 
intrinsic properties: the intrinsic absorbing column density 


(iVn), the intrinsic photon index (F), and the intrinsic X-ray 
luminosity (Lx)- Additionally, we investigate the low en¬ 
ergy X-ray spectrum of SDSS J1713+5729. The X-ray spec¬ 
tral fitting is performed using XSPEC version 12.8. lj (jArnaud 
|1996|. In all cases we account for Galactic absorption using a 
PHABS multiplicative component, with column densities fixed 
at values from |Kalberla et ak| ( |2005| l. 

4.1.1. SDSS J121839.40+470627. 7 

SDSS J1218+4706 has the strongest NuSTAR detection in 
the 8-24 keV band, with net source counts of £> 8-24 keV = 
188 for FPMA+B. The NuSTAR data are complemented by 
relatively high quality soft X-ray data, with two long XMM- 
Newton exposures (obsIDs 0203270201 and 0400560301; see 
Table [lj. Below we analyze the broad-band (0.5-24 keV) 
NuSTAR plus XMM-Newton dataset (shown in Figure |7j). The 
modeling approach taken is similar to that adopted by G14 for 
SDSS J1034+6001, the other brightest source in the NuSTAR- 
observed QS02 sample, which has comparable photon statis¬ 
tics (5*8—24 keV = 182). We group the data by a minimum of 
20 counts per bin, and use x' 2 minimisation (statistic chi 
in XSPEC) to constrain parameters. We note that using, in¬ 
stead, statist ic cst at (applying the W statistic approach; 
e.g., see Section|4.1.2|i results in essentially unchanged values 
for the key best-fit parameters (F and Aji change by less than 
0.1 and 0.1 x 10 24 cm -2 , respectively, for the models tested). 
The XMM-Newton\NuSTAR cross-normalization factor, when 
left as a free parameter, converges to slightly different values 
depending on the model being tested, but is always broadly 
consistent (given the uncertainties) with the current best cal¬ 
ibration measurements of [Madsen et al.| ( 2015| ) of ss 0.93. 
We therefore fix the cross-normalization factor to this value 
throughout. 

As shown in Figure [6] SDSS J1218+4706 has an ex¬ 
tremely flat effective photon index over the NuSTAR band, 
T 3—24 keV = —0.15 +q' 4 5. This is indicative of a spectrum 
dominated by Compton reflection, as a result of the primary 
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continuum being heavily suppressed by CT levels of photo¬ 
electric absorption (e.g., George & Fabian||1991|). Another 
important diagnostic feature of reflection is fluorescent Fe Ka 
line emission, which occurs at rest-frame 6.4 keV and be¬ 
comes incre asingly promin ent as the level of absorption in¬ 
creases (e.g.,|Risaliti|2002)>. An equivalent width threshold of 
EW Fe Ka > 1 keV is commonly used to identify CT AGNs; 
such high values are difficult to explain for less than CT 
columns (e.g., Maiolino et al. 1998| Comastri |2004] ), and sug¬ 
gest a heavily reflection-dominated or pure reflection spec¬ 
trum, where little to none of the directly transmitted AGN 
emission is visible. 

For SDSS J1218+4706, there is a clear excess of emission 
at observed frame ss 6 keV, which has previously been in¬ 
terpreted as Fe Ka line emission (J13; LaMassa et al. 2012| l. 
To model this, we fit to the > 2 keV NuSTAR plus XMM- 
Newton dataset an unobscured power law and Gaussian com¬ 
ponent, fixing the line energy at £); ne = 6.4 keV and the line 
width at (7 ii n e = 0.01 keV. We measure an observed-frame 
equivalent width of EWFe Ka = 1-7g keV using the XMM- 
Newton spectra. This value is similar to but more tightly con¬ 
strained than that published by J13, since they only use one 
of the archival XMM-Newton observations, while we use two 
here. The Fe Ka line equivalent width is above the commonly 
adopted threshold for CT AGNs (EWp e Ka > 1 keV), with a 
comparable value to that of the CT quasar SDSS J1034+6001 
(Mrk 34; G14). Freeing the Gaussian line energy parame¬ 
ter, we obtain a best-fit value of En ne = 6.40/qq 4 keV (rest 
frame), which adds further confidence that the excess emis¬ 
sion is due to Fe Ka. 

For the X-ray spectral modeling of SDSS J1218+4706, we 
first conduct a simple test to assess the nature of the AGN con¬ 
tinuum; we fit the 7-24 keV NuSTAR data with two extreme 
models, one reflection-only spectrum and one transmission- 
only spectrum. Fitting the high energy data above 7 keV al¬ 
lows a clean measurement of the AGN continuum indepen¬ 
dent of how the potentially complex lower energy emission 
is chosen to be modeled; low energy X-ray emitting pro¬ 
cesses other than the reflected or directly transmitted AGN 
continuum can dominate up to energies of « 4 keV (e.g., 
Gandhi et al.|2014||2015J>, and fluorescent line emission (e.g., 
Fe Ka) can also strongly contribute at energies up to 7 keV. 


For the reflection-only model we use PEXRAV (Magdziarz & 
Zdziarski]|1995j>, with the reflection scaling factor set to — 1 


to produce a reflection-only spectrum (i.e., no directly trans¬ 
mitted component), and set all other parameters to the default 
values. This model provides a statistically acceptable fit to 
the NuSTAR data (\ 2 /n = 11.3/12; here, n is the number 
of degrees of freedom), and the intrinsic photon index is con¬ 
strained to be T = 1.35 ± 0.46. For the transmission-only 
model we use CABS • ZWABS ■ P0W (in XSPEC formalism)]/] It 
is not possible to simultaneously constrain Nr and T in this 
case, so we fix the intrinsic photon index at T = 1.8 (a typ¬ 
ical value for AGNs detected by NuSTAR at 3-24 keV; e.g., 
Alexander et a l.|2013 1. Again, there is a statistically accept¬ 


able fit to the data Tx 2 /n = 10.5/12), for a best-fit column 


3 The model PLCABS lYaqoobj 1997) is generally a preferable trans- 
mission model to use (over CABS • ZWABS • P0W) for column densities of 
TVh > few xlO 23 cm -2 . However, in our case PLCABS is not appropri¬ 
ate, since the energy range for which the model is valid depends on source 
column density (E < 1 4.4 keV for 7 Vh < 10 24 cm -2 ; E < 10 keV for 
ATh < 5 X 10 24 cm -2 ;|Yaqoob| 1997f, which means not utilising the high 
energy NuSTAR data. 



FIG. 7.— Unfolded NuSTAR plus XMM-Newton spectrum, in EFe units, 
for SDSS J1218+4706. The data are shown in the upper panel, grouped to 
a minimum significance of 2cr per bin for visual purposes. The NuSTAR 
FPMA and FPMB data are shown in purple and orange, respectively. The 
MOS (obsID 0203270201) and PN (obsID 0400560301) data are shown in 
green and blue, respe ctively . The best-fit MYTORUS-based model (Model M; 
described in Section |4.1.1) is shown binned to match the data (solid lines, 
upper panel) and in full detail (lower panel). 

density of A H = (1.9±g;^) x 10 24 cm" 2 . 

The above tests support the empirical evidence (from r e g 
and EWpe Ka) that extremely large, CT column densities are 
required to explain the X-ray spectrum of SDSS J1218+4706. 
In the most extreme case, the source is consistent with be¬ 
ing fully reflection-dominated (no directly transmitted com¬ 
ponent), which would imply An 3> 1.5 x 10 24 cm -2 . In the 
least extreme case, the source is consistent with lying close 
to the CT threshold (An m 1.5 x 10 24 cm -2 ). However, 
the latter model assumes a transmission-only spectrum (no 
Compton reflection), which is unlikely given the large mea¬ 
sured equivalent width of Fe Ka. The reflection-only model 
tested (PEXRAV) is also limited in that the geometry (a slab 
of material) and infinite optical depth assumed are not well 
motivated for obscured AGNs. Ideally, in the CT regime, 
any absorbed continuum, reflected continuum and fluores¬ 
cent lines should be modeled in a self-consistent way, and 
assuming a well-motivated geometry. This is possi ble rising 
the physical models MYTORUS (jMurphy & Yaqoob||2009) l 
and BNTORUS (Brightman & Nan dra|201 1) , which were pro¬ 
duced using Monte Carlo simulations of X-ray radiative trans¬ 
fer through toroidal distributions of gas, with the two models 
assuming different toroidal geometries. We proceed to anal¬ 
yse the broad-band (0.5-24 keV) XMM-Newton plus NuSTAR 
spectrum of SDSS J1218+4706 using these two models. 

Our MYTORUS-based model (Model M hereafter) has the 
following form: 

Model M = PHABS x (MYTZ X P0W + MYTS + MYTL + APEC). 


Here, MYTZ reprocesses the zeroth-order transmitted contin¬ 
uum (P0W) through photoelectric absorption and the Compton 
scattering of X-ray photons out of the line-of-sight, MYTS is 
the scattered/reflected continuum produced by scattering X- 
ray photons into the line of sight, and MYTL is the fluores¬ 
cent emission line spectrum (Murphy & Yaqoob| |2009 1. We 
use MYTORUS in the simplest form possible, tying the com¬ 
mon parameters of MYTZ, MYTS and MYTL (Ah and 9 lnc ) to- 
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TABLE 4 

Best-fit Models for the NuSTAR + XMM-Newton Spectrum of 
SDSS J1218+4706 



ModelM 

ModelT 

X 2 /n 

31.9/38 

33.0/39 

r 

Z - 4 -0.3 

2-8/S.4 

N H HO 24 cm -2 ) 

zu -0.8 

9 9+1-2 
z,z -0.6 

6>tor (°) 

[60.0] 

[60.0] 

dine (°) 

63.7/*;® 

[87.0] 

feTAPEc (keV) 

0 42+ 0,20 

u -^ z -o.n 

0 2‘5 +u,uV 
u - zo -0.05 

Tf. P f-2 k eV a0 41 ergs -1 ) 

1.38 

1.65 

T^iokev CO 44 ergs -1 ) 

0.01 

0.01 

^10— 40keV a0 44 ergs -4 ) 

0.14 

0.13 

^iOkeV (!0 44 “g^ 1 ) 

0.85 

1.70 

+0-40keV a0 44 ergs -4 ) 

0.46 

0.48 


NOTE. - Best-fitting model parameters for the 0.5-24 keV spectr um of 
SDSS J1218+4706. The individual models are detailed in Section 14.1.11 
The column densities (TVh ) quoted are defined along the line-of-sight of the 
observer. 



Fig. 8 .— Intrinsic photon index (F) versus (line-of-sight) column density 
(;Y|.[) confidence contours for SDSS J1218+4706. The contours outline the 
68%, 90% and 99% confidence regions, and the best-fit value is marked by 
a black cross. We show results for two models (Model M and M odel T ; left 
and right panels). The individual models are detailed in Section |4.1.1| The 
gray shaded region indicates the parameter ranges for which Model M is not 
valid. The best-fit column densities are CT ((Vjj > 1.5 X 10 24 cm -2 ), and 
the 90% CL lower A+i limits lie just below and just above the CT threshold 
(gray dashed line) for Model M and Model T, respectively. 


gether. The intrinsic (unprocessed) photon indices and nor¬ 
malizations are tied to those of the zeroth-order continuum 
(P0W). The torus opening angle ($ t or) is fixed at 60° in the 
current version of MYTORUS. APEC is a thermal plasma com¬ 
ponent (Smith et al. |2001 ) which we use to parameterize the 
low energy excess, fixing the abundance parameter at solar. 
This component is motivated by the steep spectral slope at 
low energies (To. 5 _ 2 keV ~ 3.4, measured using an unab¬ 
sorbed power law model), which suggests contributions from 
processes such as star formation or AGN photoionization, al¬ 
though we lack the spectral detail required to distinguish be¬ 
tween these processes. The best-fit model has x 2 /n = 32/38 
(see Table [4] for the model parameters and Figure [7] for the 
model spectrum). Since T and Ah are known to be de¬ 
generate, we compute their uncertainties from \ 2 contours 
in the T-Ah plane. Contours showing the 68%, 90% and 
99% confidence regions for this parameter space are shown 
in Figure [8] These were computed with 9 nir left free to vary. 
Hereafter, the quoted uncertainties for An and T are taken 
from the 90% CL contours. The best-fit intrinsic photon in¬ 


dex and line-of-sight column density are T = 2.4//'/ and 
Ah = (2.0/q 8 ) x 10 24 cm -2 [corresponding to an equato¬ 
rial column density of AH, eq = (4.2 /q 8 ) x 10 24 cm -2 ] for 
the best-fit inclination angle of 9 mc = 63.7+2 g°- The mod¬ 
eling will not allow inclination angles of 0 UU . < 60°, since 
for these angles the observer has a direct, unobscured view of 
the central X-ray emitting source. The upper error on Ah 
is not constrained, which is in part due to the limited Ah 
range of MYTORUS (A H = 10 22 -10 25 cm -2 ). The best- 
fit model spectrum is reflection-dominated, with the MYTS 
component dominating at « 3-10 keV, and the MYTZ • P0W 
and MYTS components contributing equally to the normal¬ 
ization and spectral shape at > 10 keV. To assess whether 
the NuSTAR plus XMM-Newton spectrum is in agreement 
with being fully reflection dominated, we test two modifi¬ 
cations of Model M where the MYTZ • P0W component is re¬ 
moved and the inclination angle of the MYTS component is 
set to 0° and 90°, corresponding to face-on and edge-on re¬ 
flection. Both models provide statistically acceptable fits to 
the spectrum (x 2 /n = 29/35 and 28/35, respectively), with 
flat x 2 residuals, reasonable best-fit intrinsic photon indices 
(T = 1.6l/' G and 1.9+“, respectively) and large column den¬ 
sities for the reflecting material [An^eflector = (3.l/“ 6 ) and 
(1-5+J'g) x 10 24 cm -2 , respectively]. The broad-band X- 
ray spectrum of SDSS J1218+4706 is therefore in agreement 
with being fully reflection dominated. Since no transmission 
component is required in these models, we may infer that the 
line-of-sight column density is consistent with having a value 
of A h 3> 1.5 x 10 24 cm~ 2 . 

Our BNTORUS-based model (Model T hereafter) has the 
following form: 

Model T = PHABS X (BNTORUS + APEC). 


In the BNTORUS model. Ah is defined along the line of sight, 
and is independent of 0\ nc . Initially, we fix the inclination at 
the maximum value of 9 lnc = 87°, corresponding to an edge- 
on view of the torus. Since the opening angle for Model T is 
poorly constrained when left as a free parameter (9 to r < 72°), 
we fix it to 60°. The best-fit model has \' 2 /n = 33/39 (the 
model parameters are listed in Table [4] and the T-Ah con¬ 
tours are shown in Figure [8j. Ah is well constrained at the 
90% CL, with a best-fit value of (2.2+g'g) x 10 24 cm -2 , 
and the intrinsic photon index has a relatively high value of 
T = 2.8 +q 4 . The upper error on T is not constrained due to 
the parameter limits of the BNTORUS model. Fixing the in¬ 
trinsic photon index at a more reasonable value of T = 2.3, 
which is consistent with the x 2 contours and is at the up¬ 
per end of the range typically observed for unobscured AGNs 
(e.g., Mateos et al.|2010||Scott et al. 201 l| l, results in a higher 
column density of Ah = (3.6/J/S) x 10 21 cm -2 and a re¬ 
duced x 2 value close to unity (x 2 /n = 39/40). If the in¬ 
trinsic photon index is fixed at T = 1.8, an extremely high 
column density of Ah > 5.1 x 10 24 cm -2 is required. We 
note that the modeling (with T left free) allows a large range 
of inclination angles ((/ tK: > 63°), and re-modeling with 0\ nc 
fixed at a lower value of 65° results in a similarly good fit 
(X 2 / n = 38/39) with no significant change in Ah but a flat¬ 
ter photon index of T = 2.5/+ij. Furthermore, the statistical 
quality of the fit and the best-fit parameters are relatively un¬ 
changed when (Lor is left as a free parameter. 

To summarize, CT line-of-sight column densities are pre¬ 
ferred for all of the models tested for SDSS J1218+4706. 
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Fig. 9.— Unfolded NuSTAR (purple and orange for FPMA and FPMB, re¬ 
spectively) plus Chandra (blue) X-ray spectrum for SDSS J1243-0232. The 
best-fit absorbed power law (ZWABS • POW) model is shown. The panel layout, 
units and data binning follow that of Figure^] 


The broad-band X-ray spectrum shows evidence for having 
a dominant contribution from Compton reflection, with the 
primary continuum being heavily suppressed due to photo¬ 
electric absorption. This is in agreement with the expec¬ 
tations from the observation of strong fluorescent Fe Ka 
line emission (EWp e Ka ~ 1.7 keV). The lowest limit al¬ 
lowed by the modeling for the line-of-sight column density is 
JVh > 1.2 x 10 24 cm -2 , and there is no constraint at the upper 
end. The JVh, £x and EWp e Ka constraints and data quality 
for SDSS J1218+4706 (z = 0.094) are remarkably similar 
to those for the other low redshift QS02 strongly detected 
by NuSTAR , SDSS J1034+6001 (z = 0.051; also known as 
Mrk 34), which was identified by G14 as a bona fide CT AGN. 
More complex models are possible (such as a clumpy torus; 
e.g., Bauer et al.|2014[ ), but testing these is beyond the X-ray 
data quality. 


4.1.2. SDSS J124337.34-023200.2 

SDSS J1243-0232 is the third brightest NuSTAR detec¬ 
tion in the SDSS-selected candidate CTQS02 sample, after 
SDSS J1218+4706 (Section [4~TTT| ) and SDSS J1034+6001 
(G14), but still has relatively low photon counts: Ss-24 keV ~ 
90 and S3-8 keV ~ 34 with NuSTAR, and So. 5-8 keV ~ 9 
with Chandra. This emphasizes the challenge involved in 
studying these inherently faint X-ray sources. Due to the low 
photon statistics, we use statistic cstat in XSPEC, which 
is more appropriate t han statistic chi in the case of Pois¬ 
son distributed data ( |Nousek & Shue[[l989[ ). In the case of 
un modeled background spectra, cstat applies the W statis¬ 
tic fWachter et al. 1979) J While the W statistic is intended 
for unbinned data, bins containing zero counts can lead to er¬ 
roneous results]^] so we group the Chandra and NuSTAR data 
by a minimum of _1_ count and 3 counts per bin, respectively 
(e.g., |Wik et al.||2014 i. We fix the Chandra.NuSTAR cross¬ 
normalization factor at 1.0, consistent with the value obtained 
when the cross-normalization factor is left as a free parameter 
in the modeling. 


4 See also http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/wstat.ps 

5 See https://heasarc.gsfc.nasa.gov/xanadu/xspec/ 
manual/XSappendixStatistics.html 


The NuSTAR spectrum of SDSS J1243-0232 has a flat ef¬ 
fective photon index of T3_24 keV = 0.66±0.50, indicative of 
heavy absorption. Fitting the broad-band (0.5-24 keV) NuS¬ 
TAR plus Chandra spectrum with a simple absorbed power 
law (ZWABS • POW) model, we obtain JV H ~ 1.6 x 10 24 cm -2 
and T « 3. This intrinsic photon index is discrepant with the 
expected range for AGNs, and the parameter is poorly con¬ 
strained. We therefore fix the parameter to T = 1.8 (typical 


value in the 3-2 4 keV energy band for AGNs; e.g., Alexan-| 
der et al.]|2013 1. The best-fitting model has x 2 = 101 and 


a C-statistic value of C = 123, for n = 130. The un¬ 
folded spectrum and best-fitting model are shown in Figure 
[ 5 ] The column density, JVh = (0.90+gp) x 10 24 cm -2 , is 
close to CT. The intrinsic luminosities in the low and high 
energy X-ray bands are C^iokeV = 9.6 x 10 44 erg s ” 1 

and £' 10-40 keV = 9.7 x 10 44 erg s _1 , respectively. The 
higher quality NuSTAR data dominate the fit, with similar re¬ 
sults [JV H = (0.97^Q 3 g) x 10 24 cm" 2 ] being obtained when 
the Chandra data are excluded. We note that cstat may also 
be used to model the unbinned, gross (i.e., combined source 
plus background) spectrum, in which case the Cash statistic 
(C statistic; |Cash|1979 1 is applied. Characterizing the back¬ 
ground spectra using double power law models (POW + POW in 
XSPEC), and including these as fixed components in the spec¬ 
tral modeling of the NuSTAR data, this C statistic approach 
yields very similar results to the W statistic approach, with 
JV H = (0.97+°;p) x 10 24 cm" 2 . 

Given the extremely flat effective photon index measured 
for this source, it is reasonable to test whether the spectrum 
is i n agree ment with a pure reflection continuum. As in Sec- 
tion |4. 1. l] we use PEXRAV with the reflection scaling factor set 
to —1 to produce a reflection-only spectrum. The model pro¬ 
duces a similarly good fit to the data as for the absorbed power 
law model above, with y 2 = 117 and C = 120, for n = 130. 
We infer that the line-of-sight column density is consistent 
with being CT, with JVh 1.5 x 10 24 cm -2 . Unlike for the 
absorbed power law model, the intrinsic photon index is well 
constrained by the reflection-only model, with T = 1.7 ± 0.3. 
To summarize, the NuSTAR data unambiguously reveal heavy 
absorption in this QS02, with a column density lower limit of 
JVh > 0.6 x 10 24 cm -2 and no constraint at the high, CT ab¬ 
sorption end. Higher quality X-ray data than those currently 
available, especially at < 10 keV, are required to reliably dis¬ 
tinguish between less than CT, and reflection-dominated CT 
models. For instance, the current data are unable to provide 
informative constraints on Fe Ka line emission (see the Ap¬ 
pendix). 


4.1.3. SDSS J171350.32+572954.9 


For SDSS J1713+5729 there are too few NuSTAR counts 
for broad-band X-ray spectral modeling (see Table [2]). Here 
we investigate the low energy (< 10 keV) spectrum observed 
with XMM-Newton. The object appears to have an extremely 
steep spectrum at low energies, with PN (MOS) source counts 
of < 2 (< 5) at 2-10 keV and 12±| (18±|) at 0.5-2 keV, 

implying a photon index of T = 3.5+o'g in the 0.5-10 keV 
energy band; J13 measure a slightly flatter, but consistent 
(within the uncertainties), value of F = 2.5 ± 0.4. The steep 
spectral slope is not typical of an AGN, and would be incon¬ 
sistent with the NuSTAR detection if produced as a result of 
direct AGN emission. To test whether the soft X-ray emis¬ 
sion could be powered by star formation, we compare the 0.5- 
8 keV luminosity, £0.5-8 keV = 1-4 x 10 42 erg s _1 , with the 
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Fig. 10.— NuSTAR band ratio (BRn u ) versus z. The circles, the diamond 
and the triangle indicate the candidate CTQS02s presented in this work, L14 
and G14, respectively, which are detected at 8-24 keV; 90% CL error bars 
and limits are shown. For comparison, gray s quares show the first10 AGNs 
detected in the NuSTAR serendipitous survey (Alexander et al.|2013); 68% 
CL error bars are shown. The gray pentagon marks a Cl AGIN identified 
with NuSTAR in the COSMOS field (F Civano et al., submitted); a 90% CL 
lower limit is shown. The tracks show model predictions for BRnu for four 
absorbing column densities in the range Nn = (0.1-1.5) x 10 24 cm -2 . The 
BR Nu constraints for SDSS J1218+4706, 1243-0232 and 1713+5729 are 
higher than that of the confirmed CTQS02 SDSS J1034+6001, and suggest 
large absorbing columns. 


far-infrared (FIR) luminosity, Lfir < 4.0 x 1 0 44 e rg s 1 , 
measured using IRAS fluxes following Lonsdale Persson & 


Helou (1987). The relatively high soft X-ray:FIR luminos 
ity ratio of L 0 5 _ 8 kev/^FiR > 0.0035, which is a conser¬ 
vative lower limit due to the poorly constrained IRAS 100//m 
flux, rules out star formation as the driver of the soft X-ray 
emission (e.g., see Figure 8 of Alexander et al.| |2005) . We 
deduce that the soft X-rays detected with XMM-Newton are 
indirectly powered by the AGN (e.g., via photoionization or 
scattered AGN emission), and NuSTAR may have provided 
the first identification of the directly transmitted (or reflected) 
AGN continuum of this QS02. 


4.2. X-ray Spectral Analysis: Band Ratios 

X-ray band ratios provide a basic description of the X-ray 
spectrum, and are useful when there are insufficient counts 
for detailed spectral modeling. We define the NuSTAR band 
ratio (BRn u ) as the ratio of net source counts in the hard- 
band to those in the soft-band, S’s-24 kev/^-s keV- Figure 
TO] shows BRnu against redshift (z) for the five (of the to¬ 
tal nine) NuSTAR-observed candidate CTQS02s which are 
detected at 8-24 keV, including the three presented in this 
paper (SDSS J1218+4706, 1243-0232 and 1713+5729) and 
the two presented in L14 and G14 (SDSS J0011+0056 and 
1034+6001, respectively). The tracks show the expected evo¬ 
lution of BRnu with z for four different fixed column den¬ 
sities (Nh), computed using a MYTORUS model with an 
intrinsic photon index of T = 1.8. We compare the mea¬ 
sured BRnu values for the candidate CTQS02s with these 
tracks in order to infer A r H . We note that producing the 
tracks with, instead, a simple ZWABS ■ P0W model results in 
higher Nh values for the same liR\j u . The NuSTA /^-detected 
candidate CTQS02s, in general, have high band ratios com¬ 
pared to AGNs detected in the NuSTAR extragalactic surveys 


(squares in Figure ]T0|. In all cases the BRn u values suggest 
Nh > 10 23 cm -2 . 

For SDSS J1713+5729, a NuSTAR-detected object with too 
few counts for broad-band s pectra l modeling of the direct 
AGN continuum (see Section ft. 1.3| ), the lower limit in BRn u 
suggests heavy absorption with Nh > 5 x 10 23 cm -2 . Our 
most direct measurement for the intrinsic X-ray luminosity 
of this QS02 comes from using this Nh constraint. Taking 
the observed 10-40 keV luminosity constraint from Table [3| 
and assuming that the X-ray spectrum is an absorbed power 
law with r = 1.8, the lower limits obtained are L™_ w keV 
> 4.6 x 10 42 erg s -1 and L 1 1 I J ) _ 40 keV > 5.3 x 10 42 erg s -1 . As 
an alternative to the I illy,, approach, Nh can be constrained 
using the NuSTAR/XMM-Newton band ratio (following L14). 
However, in this case the constraint (Nh > 2 x 10 23 cm -2 ) is 
less stringent than that from BRn u , due to the comparatively 
poor quality of the available XMM-Newton data. 

The Nh estimates made from BRy u using Figure [TO] are 
relatively crude, since the individual X-ray spectra mayhave 
additional spectral complexities (e.g., line emission around 
« 6.4 keV, a scattered power law, or a complex absorber ge¬ 
ometry) not incorporated in our model predictions. To illus¬ 
trate this, for the two sources with comparatively high qual¬ 
ity NuSTAR spectra (SDSS J1034+6001 and 1218+4706), the 


less than CT column densities inferred from the I Uly,, anal¬ 
ysis (Nh < 5 x 10 23 cm -2 and < 10 24 cm -2 , respectively) 
are an underestimate of the column densities determined from 
X-ray s pectra l fitting (Nh > 1.5 x 10 24 cm -2 ; see G14 and 
Section |4.1.1 of this paper, respectively). Similarly, using the 
NuSTAR results for three CT reflection-dominated Seyfert 2s, 
Balokovic et al.](J2014|) demonstrate that the above RRy u ap¬ 


proach underestimates Nh for reflection-dominated AGNs. 
Nevertheless, BRn u provides first-order Nh constraints for 
weakly detected sources. 


4.3. Indirect Constraints on X-ray Absorption 

It is well-established that there is a tight relation between 
the mid-IR and intrinsic X-ray luminosities of AGNs (e.g.. 
Lutz et al.| 12004] |Fiore et al.||2009l |Gandhi et all [2009] 

Lanzuisi et al.||2009| [Mateos et al.|2015[|Stem|2015) . Mid- 


IR emission can therefore provide an indirect estimate of 
the intrinsic AGN power, especially useful when heavy ab¬ 
sorption in the X-rays makes this information challenging 


to obtain (e.g.. 

Vignali et al. 

2010; Alexander et al. 2008 

LaMassa et al. 

2009 5011 G 

oulding et al. 2011; 

Lanzuisi| 

et al. 2015b). 

following the 

approach used for ot 

her NuS- 


TAR studies of faint, obscured AGNs (L14; |Stern et al.|2014 1, 
in Figure m we compare the observed X-ray:mid-IR lumi¬ 
nosity ratios with intrinsic ratios for unobscured AGNs and 
those corresponding to X-ray absorption due to dense obscur¬ 
ing material (Nh = 10 24 cm -2 ), for both the low (2-10 keV) 
and high (10-40 keV) energy X-ray regimes. We show the 
full sample of nine NuSIA /^-observed SDSS-selected candi¬ 
date CTQS02s, including the five presented in this work, the 
three from L14 and the one in G14. The X-ray luminosities 
(Ly hs ) are observed values (i.e., uncorrected for absorption), 
and the 6 pm luminosities ( Lq ^m- in units) are intrinsic 
values (i.e., corrected for dust extinction occuring in the sys¬ 
tem) for the AGN determined through SED modeling (Sec- 
tion |3.3[ ), and both correspond to the values provided in Table 
[3] We note that for a large fraction of CT AGNs, potentially 
« 50% in the case of local CT AGNs, we e xpect significant 
absorption in the mid-IR (e.g., |Bauer et al.||2010] |Goulding 
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Fig. 11.— Observed (i.e., uncorrected for absorption) X-ray luminosity for the rest-frame 2-10 keV and 10-40 keV bands (left and right hand panels, 

respectively) versus rest-trame 6 jin l luminosity (in uL u units). 1 he circles indicate the A/u$lAR- observed SDSS-sefected candidate CTQSU2s presented in this 

work, L14, and G14 (z = 0.05-0.49); colored circles mark the NuSTAR-detected sources. The X-ray luminosities for the candidate C lOS0"2s are taken from 

best-fitting spectral modefs"where possible. Otherwise, they have been determined from photometry, assuming an unabsorbed power law model with f ’ = 0.3 

(as described in Section|3^. For the three L14 objects (SDSS J0011+0056, 01356+0032 and 1157+6003), the values have been adiusted tor consistency with 

this work. Other AwAfAR-observed objects are shown, including: A luSTAK extragalactic survey AGJNs d 1 -!-’ symbols; z = 0.02-2.92;(Alexander et al. 2013), 

three cTT 1 Seyfert 2 AGNs (‘X' symbols; 2 : ps 0.01; Balokovic et al. 2014 1 , a heavily obscured quasar identified in the FCDFS field (star; 2 : ps 2; De Moro 

et al. 2014 1 , three luminous and heavily obscured VkAS’E-selected AGNs (diamonds; z p^^J; Stern et al.|2U14>, and a CT ACjft identified in the COSMOS held 

(pentagon; z = U.U44; F Uivano et al., submitted). For the latter object, we show an upper limit in Lq uui , since we have assumed that the mid-fR emission is 

AGN-dominated. Additionally, tor tne 2-10 keV band we compare witn sources studied at < 10 KeVwitn Chandra or XMM-Newton: tne candidate Cl quasars 

presented in Alexander et al. 12008, squares; z p^ 2), and a candidate heavily Cl AGJ^Tidentified in the COSMOS field (triangle; z = 0.35;|Lanzuisi et al. 

2015b 1 . For four of the NuSTAR^observed candidate CTQS02s, vertical lines indicate the intrinsic (i.e., corrected tor absorption) X-ray luminosities obtained 

from A-ray spectral analyses. We compare afl of the data with two intrinsic relations tor the 2-10 keV band (solid black lines), those ot|Fiore et al.|(20U9| and 

Gandhi et al. (2009i. Following L14 and Stern et al. (2014 1 , the relations have been extrapolated to 10-40 keV assuming I' = and the dashed lines sfiow 

the effect of absorption by TVh = 10 Z4 cm~ z gas. The different X-ray:mid-IR ratios for the Fiore et al. (2009) and Gandhi et al. ( 2009) relations means that the 

former provicfeTa more conservative estimate of the CT absorption threshold. The majority of the A/uSTAR-observed candidate CTGSOzs have low X-ray:mid~lR 

ratios, suggestingTj 1 levels of photoelectric absorption. 


et al.] |2012| ). We have partially addressed this through dust 
corr ections which are included in the SED modeling (Section 
|3.3| >. These corrections are small, however, with the lumi¬ 
nosities changing by factors ranging from 1.03 to 1.46 (with 
a median of 1.17). For the four candidate CTQS02s with 
constrained intrinsic X-ray luminosities (L 1 ^), we plot the 
L x* values obtained from X-ray spectral analyses (see L14, 


G14, and Sections 4.1.1 and |4. 1.2] of this work). We conser¬ 
vatively adopt intrinsic X-ray luminosities from the models 
with lower best-fit column densities (e.g.. Model M in the case 
of SDSS J1218+4706 and the absorbed power law model in 
the case of SDSS J1243-0232). 

The two intrinsic re lations ut ilized for c omparison are those 
of |Fiore et al.| ( |2009| > and Gandhi et al.| ( |2009] >, which were 
both c omputed at 2-10 keV. In the case of the |Gandhi et al. 
( |2009[ ) relation, we adjust the 12 //m (the mid-IR wavelength 
at which the relation was computed) vL v luminosities down 


wards by 7% to obtain 6 /im luminosities, based on the Assef 
et al. (2010) AGN template. The two relations predict slightly 


different X-ray:mid-IR ratios at low luminosities and diverge 
further towards higher luminosities, which is partly due to the 
different luminosity ranges over which the two relations were 
calibrated, but also reflects the uncertainty in such relations. 
Comparison to both allows us to account for systematic ef¬ 
fects in the derivation of these relations. We extrapolate the 
relations to the 10-40 keV band assuming T = 1.8 (typical 
value for AGNs; e.g., |Alexander et al .|2013j >. An advantage 
of using 10-40 keV X-ray luminosities (T°q 1 40 koV ), as op¬ 
posed to 2-10 keV luminosities (T^io key)’ t ^ lat contami¬ 


nation from processes other than AGN continuum emission is 
negligible in this high-energy band. However, the suppression 
of the X-ray emission by absorbing gas is less dramatic in the 
10-40 keV band, as demonstrated by the relative normaliza¬ 
tion of the Nh = 10 24 cm -2 lines in the left and right hand 
panels of Figure [IT] which were computed assuming a MY- 
Torus model witnT = 1.8 and (),,}, K = 70° (following F14). 
Absorption by TVh = 10 24 cm -2 gas results in a suppression 
of the X-ray emission by factors of « 20 and « 2 in the 2- 
10 keV and 10-40 keV bands, respectively. We note that for 
the four candidate CTQS02s with L'y' values constrained us¬ 
ing X-ray spectral analyses, the intrinsic luminosities agree 
mo re closely with the|Gandhi et al.| ( |2009| l relation than with 
the |Fiore et al.| ( |2009l > relation. 

In general, the overall sample of candidate CTQS02s have 
extremely low 2-10 keV:mid-IR ratios, with the observed 2- 
10 keV luminosities a factor of > 20 lower than the intrinsic 
relations, suggesting CT absorption. This was already appar¬ 
ent from 2-10 keV luminosities published in the literature, but 
here we have demonstrated the 2-10 keV suppression using 
our own soft X-ray analysis. A similar conclusion is reached 
in the high-energy 10-40 keV band, where six out of nine 
of the objects have X-ray luminosities a factor of > 2 lower 
than the intrinsic relations, consistent with CT obscuration. 
Our sample of SDSS-selected candidate CTQS02s lies below 
the majority of the AGNs detected in the NuSTAR extragalac- 
tic surveys (Alexander et al. 2013[ >, including a heavily ob¬ 
scured quasaFHetecteHTn^CDFSTNuS'lAR J033202-2746.8; 
z ~ 2; |Del Moro et al.|20141 i. 

Of the five new objects presented in this work, there is one, 
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SDSS J1243-0232, which does not appear compatible with 
CT absorption based on this indirect analysis. For this object, 
the low JVh implied by the relatively high X-ray:mid-IR ratios 
is incongmous wit h the d irect constraints from X-ray spectral 
modeling (Section 4.1.2 1 , which suggest Ah > 10 24 cm -2 . 
A similar case where the Ah values inferred from X-ray spec¬ 
tral modeling and the X-ray:mid-IR ratio do not agree is that 
of NuSTAR J033202-2746.8 (star symbol in Figure fTT] |DeI 
Moroet al.|2014l. Despite the large co lumn density me asured 


- 2 . 


Del Moro et al. 2014}, 


for this source (7Vh « 6 x 10 23 cm 
it lies high with respect to the relations, which may in part 
be due to its significant Compton reflection component. It is 
possible that a strong reflection component also contributes to 
the high X-ray:mid-IR ratio observed for SDSS J1243-0232, 
especially given that a pure reflection spectrum well describes 
the data (see Section [4.1.2[ ). 

Of the NuSTAR targets detected at high energies (> 
10 keV), SDSS J1713+5729 has the most extreme 10- 
40 keV:mid-IR ratio, with a L°pl 40 kcV value suppressed 
by a factor of ss 35 with respect to the intrinsic relations 
(on average). The fact that the source lies even lower 
than the CTQS02 SDSS J1034+6001 (G14) may be due 
to some combination of a heavily CT absorbing column 
(iVf [ 3g> 10 24 cm -2 ) and a less prominent reflection com¬ 
ponent. For the non detections, SDSS J0758+3923 and 
SDSS J0840+3838, the L°q1 40 keV upper limits suggest that 
if the X-ray faintness is due to absorption, these sources are 
likely CT (for SDSS J0840+3838 this only applies for the 
Gandhi et al. [2009| relation). While heavy absorption seems 
the mosflikely explanation for the X-ray faintness of these 
non detections, we do not have broad-band X-ray spectral 
constraints and therefore cannot rule out the possibility of in¬ 
trinsic X-ray weakness (e.g., Gallagher et al.|200l] Wu et al. 


2011 Luo et al. 2014| Teng et al.||2014[ ). However, intrinsic 

X-ray weakness is a phenomenon observed for type 1 sources 
where there is an unobscured view of the central nucleus, un¬ 
like for our QS02s. 

5. DISCUSSION 

In the following sections, we discuss the possible impli¬ 
cations of the extremely high column densities and corre¬ 
sponding intrinsic luminosities measured for the NuSTAR- 
detected heavily obscured QS02s presented in this pa¬ 
per (SDSS J1218+4706, 1243-0232 and 1713+5729), L14 
(SDSS J0011+0056), and G14 (SDSS J1034+6001), in the 
context of the overall quasar population. 

5.1. Heavy Absorption and Powerful X-ray Luminosities 


Figure 12 shows N h versus intrinsic (i.e., absorption- 
corrected) X-ray luminosity for all SDSS-selected QS02s that 
have been studied at low energies (< 10 keV) with Chan¬ 
dra and XMM-Newton, and have direct constraints from X-ray 
spectral analyses. The intrinsic X-ray luminosities shown are 
for the rest-frame 2-10 keV band (L]f_ 10 koV ), and are here¬ 
after referred to as Lx - The data are compiled from J13 and 
|LaMassa et al.| ( |2014| >. Since these two studies have differ¬ 
ent approaches, with the former limiting the spectral analysis 
to absorbed power law models and the latter using physically 
motivated models, we adopt the LaMassa et al.] ( |2014| val¬ 
ues where multiple measurements exist. Overlaid are the five 
sources which have 8-24 keV detections with NuSTAR , for 
which it is therefore possible to remeasure JVh and Lx with 
the addition of the high-energy (> 10 keV) data. In each case, 
there is a range of column densities consistent with the data. 



Fig. 12.— Intrinsic (i.e., absorption-corrected) rest-frame 2-10 keV lumi¬ 
nosity (IAf j q keV , or Lx) versus JVh for SDSS-selected QS02s, as mea¬ 
sured from X-ray spectral analyses. The grey open circles show the con¬ 
straints in the literature, all directly measured from X-ray spectral fitting at 
soft X-ray energies (< 10 keV; J13;|LaMassa et al.|2014|. Sources with ev¬ 
idence for being CT, primarily based on the indirect X-ray: [O III] luminosity 
ratio diagnostic, are marked with a “C”. The colored circles show our con¬ 
straints for the five NuSTAR-observed candidate CTQS02s detected at high 
energies (> 10 keV), from the broad-band NuSTAR plus soft X-ray spectral 
analyses presented in this study, L14, and G14. The colored lines indicate 
the significant increase in both Lx and N h for these five objects between 
the soft X-ray constraints in the literature and the NuSTAR plus soft X-ray 
studies. We note that for SDSS J1243-0232 the increase in Lx shown (blue 
line) may be an overestimate]^ 

“For SDSS J1243-0232, we measure a significantly higher Chandra flux 
(by roughly a factor of eight) than J13 using the same data. While there is 
not a clear reason for this discrepancy, we note that our measured Chandra 
3-8 keV flux agrees well with the NuSTAR flux for the same energy band (see 
Table[3}. 


To be conservative, we adopt measured values at the lower 
end of these ranges: e.g., for SDSS J1218+4706 we adopt 
the Model M results (Ah = 2.0 x 10 24 cm -2 ; Section 4.1,1 1 
and for SDSS J1243-0232 we adopt the abso rbed p ower law 
model results (Ah = 9 x 10 23 cm -2 ; Section |4~L2 1 . The im¬ 
provements made with NuSTAR are illustrated by the colored 
lines, which connect the literature constraints prior to NuSTAR 
and the broad-band, NuSTAR plus soft X-ray constraints. 

Our Lx and Ah measurements for these five objects are 
significantly higher than the constraints in the literature from 
spectral modeling of the soft X-ray (Chandra or XMM- 
Newton) data alone. For the fainter quasars which have 
net Chandra (0.5-8 keV) or XMM-Newton PN (0.5-10 keV) 
source counts of 5 so f t < 15 (SDSS J0011+0056, 1243-0232 
and 1713+5729) the soft X-ray constraints underpredict Ah 
by factors of /cjv h ~ 30-1600, while for the brighter sources 
with S s 0 ft > 50 (SDSS J1034+6001 and 1218+4706) A H is 
underpredicted by factors of fcjv H « 2.5-5. In general, the 
intrinsic X-ray luminosities (Lx) measured are ss 1-2 orders 
of magnitude higher with the addition of NuSTAR data, which 
is largely due to the increased absorption correction. These 
results have implications for X-ray studies of AGNs at z < 1 
that lack sensitive high-energy (> 10 keV) coverage. For ex¬ 
ample, on the basis of our results we infer that X-ray data at 
< 10 keV may not reliably identify heavily obscured to CT 
(Ah > 5 x 10 23 cm -2 ) AGNs if the photon counts are low, 
and the intrinsic luminosities will be underestimated. A simi¬ 
lar conclusion was reached by [Wilkes et al.|([2013|l, who used 
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Fig. 13.— The iVjj distribution of SDSS-selected QS02s at z < 0.5, con¬ 
structed using only direct constraints from X-ray spectral analyses. The five 
NuSTAR-observed objects with enough > 10 keV counts for X-ray spectral 
analyses are marked by the hatched regions. Left panel: a measurement of 
the A'i i distribution from existing soft X-ray (< 10 keV) Chandra and XMM- 
Newton cons traints. CT candidates, identified as such in J13 and LaMass a 
|et al.|(2014| primarily using the indirect X-ray: [O III] luminosity ratio di¬ 
agnostic, are labeled as “CT”. Sources with upper/lower limits in Aj | are 
marked with arrows. Right panel: Prediction for the true Nn distribution 
of SDSS-QS02s, based on the results of this broad-band NuSTAR plus soft 
X-ray study; see Scction [l.2| 

Chandra and multiwavelength data to investigate the intrinsic 
X-ray properties of quasars selected at low radio frequencies. 

The intrinsic X-ray luminosities of our objects (close to 
Lx = 10 44 erg s -1 , which ro ughly agr ees with the Lx,* value 
for unobscured AGNs; e.g., Hasinger et al. 2005) makes them 
important for population synthesis models of the CXB, since 
z < 1.5 AGNs around this lumin osity produce most of the 
CXB at its high energy peak (e.g., Treister & Urry 2005 )£] It 
is thus useful to consider the Ah distribution and CT fraction 
for this class of optically selected QS02s. 

5.2. The Ah Distribution 

In the left panel of Figure [13] we show the observed Ah 
distribution for SDSS-selected QS02s that are detected with 
Chandra and XMM-Newton, and have direct constraints at 
< 10 keV from X-ray spectral fitting (J13; LaMassa et ah 
2014). The 39 objects included have z < 0.5 and L[ q ~ 
> 2.5 x 10 s Lq, and should therefore be broadly represen¬ 
tative of the overall optica lly s elected QS02 population (for 
further details, see Section [T2) . The exclusion of QS02s un¬ 
detected by Chandra and XMM-Newton has a negligible im¬ 
pact since, for the adopted z and L [q m i ranges, there are only 
three such objects. On the basis of these data, the column den 
sity distribution is relatively flat at Ah = 10 21 -10 24 
and there is only one object above Ah = 10 24 cm -2 . The 
absorber for this object (SDSS J0939+3553) appears different 
in nature to those presented in this paper, possibly taking the 
rare form of a geometrically thin toroidal ring (LaMassa et al. 
|20l4) . 

In this work, we have demonstrated that soft X-ray ( Chan¬ 
dra and XMM-Newton ) studies can underpredict the Ah and 
Lx values of quasars with evidence for CT absorption based 
on m ultiwavelength diagnostics (CT candidates; see Section 
ED and Figure |72). The severity of the Nn and Lx under¬ 
predictions is related to the observed soft X-ray source pho- 

6 While the A/w5T4.K-detected obj ects a ll s atisfy the classical optic al qu asar 
luminosity definition (see Sections ED and \2.2\ , based on Figure |T2] they 
are just below the standard ‘X-ray quasar’ luminosity threshold {Lx > 
10 44 erg s -1 ), although SDSS J1034+6001. 1218+4706 and 1243-0232 
are consistent with lying above the threshold for some of the X-ray spectral 
model solutions. 


2 


ton counts (SsoftX with the faintly detected sources suffer¬ 
ing larger underpredictions than the more strongly detected 
sources. To understand the consequences of this for the 
true Ah distribution of QS02s, our result for the NuSTAR- 
detected objects can be extrapolated to the remaining CT can¬ 
didates in Figure 13 which were identified as such primar¬ 


ily based on the X-ray:[O III] luminosity ratio (J13; LaMassa 
et al. |2014). This extrapolation relies on assuming that the 


AwS7AR-detected subsample of five objects are representa¬ 
tive of the remaining subsample of 19 CT candidates in terms 
of their absorption properties. This is a reasonable assump¬ 
tion; the Lx bs /Lg^m distributions of the two subsamples are 
in agreement (KS test: p = 0.70), using the X-ray luminosi¬ 
ties from J13 (except for SDSS J1243-0232, for which we use 
our measured luminosity; see footnote [a) and estimating the 
6 /mi luminosities from an interpolation between the WISE 
photometric bands. 

To make a prediction for the true Ah distribution of op¬ 
tically selected QS02s, we apply an Ah correction factor 
(fc/v H ) to each of the 19 CT candidates in Figure [13] not 
observed/detected with NuSTA R , in formed by our NuSTAR- 
measured k]y H values (Section |5.1 b. For sources with low 
(S so ft < 33) and high (S so{t > 33) soft X-ray source 
counts (using PN counts only in the case of XMM-Newton 
data) we draw correction factors at random from flat dis¬ 
tributions between 1.5 < log(/c at h ) < 3.2 and between 
0.4 < log(fcj\r H ) < 0.7, respectively. In determining these 
correction factors we assumed column densities which are at 
the lowe r en d of the range that is consistent with the data 
(Section |5.1[ ): for the three most strongly detected sources 
(SDSS JHB4+6001, 1218+4706, 1243-0232), the lowest 
best-fit Ah values of (0.9-2.0) x 10 24 cm -2 are adopted, 
although the sources are consistent with having much larger 
columns (Ah > 5 x 10 24 cm -2 ); and we assume the Ah 
lower limit for SDSS J1713+5729 (Ah = 5 x 10 23 cm -2 ). 
As such, the Ah distribution prediction below may provide 
a lower limit on the CT fraction. However, this discussion 
is ultimately limited by the small number of sources detected 
above 10 keV with NuSTAR. 

The predicted Ah distribution (averaged over many itera¬ 
tions) is shown in the right hand panel of Figure [13] This 
“NuSTA ^-informed” Ah distribution for optically selected 
QS02s is strongly skewed towards high columns of Ah > 


10 23 cm 2 . Our predicted CT fraction (/ct)> defined here 

~|24 


as the ratio of the number of objects with Ah > 10" cm 
to those with Ah > 10 22 cm -2 , is /ct = 361)2 where 
the errors represent binomial uncertainties only. The full un¬ 
certainties are likely to be larger; considering extreme k\r u 
distributions, where the values assumed are all set equal 
to either the highest or lowest values of the ranges measured 
with NuSTAR, the uncertainties on /ct may be larger by a 
factor of ss 2. 

The CT fraction is an important parameter in population 
synthesis models of the CXB. In many such models, /ct is 
treated as a fixed, global quantity; the Treister et al.; ( |2009) 
model assum es a relatively low value of 15%, w hile others 
assume 50% ( Gilli et aL]|2007 Ueda et ak]|2014 the quoted 
fractions have been adjusted from the original published val¬ 
ues to our adopted definition of /ct)- It is possible to esti¬ 
mate /ct using this class of CXB synthesis model, although 
meaningful constraints are challenging to obtain due to de¬ 
generacies with other parameters (e.g., Akylas et al. 2012]). 
Fixing the Compton-reflection strength parameter, Ueda et al. 
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(2014) constrain /ct = 33-62%, which is compatible with 
our result. In other CXB synthesis models, the CT fraction 
is dependent on physical properties of the AGN population; 
according to the |Draper & Ballantyne| ( |2010] l model, high CT 
fractions are associated (beyond the local Universe) specif¬ 
ically with black holes accreting at a large fraction of their 
Eddington rate, in broad consistency with our findings. 

With the (Yu distribution in Figure [13] we have attempted 
to provide a prediction using only directly measured column 
densities since analysis of the X-ray spectrum should pro¬ 
vide the “purest” measurement of the line-of-sight column 
density, without the need to make assumptions in comparing 
emission across very different wavelength regimes (i.e., us¬ 
ing indirect absorption diagnostics such as the X-ray:mid-IR, 
X-ray:[O III] or X-ray: [Ne V] luminosity ratios). However, 
it is worthwhile considering an extreme scenario in which 
all of the candidate CTQS02s in Figure [13] (labelled “CT”) 
are truly CT; i.e., in which the indirect absorption diagnos¬ 
tics are assumed to be accurate. Applying this assumption, 
the predicted CT fraction is /ct = 651 [3 %• For compari¬ 


son, [Wgnairerar]( 20T0]) make similar assumptions using the 
X-ray:[0 III] and X-ray:mid-IR luminosity ratios for a com¬ 
plete sample of 25 SDSS-selected QS02s at z ~ 0.5, and 
determine /ct ~ 50%. Additionally, Vignali et al. (2014 
utilize the X-ray:[Ne V] ratio for a sample of z ps 0.8 type 
AGNs and find /ct ~ 40%. In the case of Seyfert 2s in the 
local Universe, A/j distributions have been constructed for 
optically selected samples using indirect absorption diagnos¬ 
tics (primarily the X-ray:[O III] ratio), predicting a fraction of 
/ct 50% for this lower luminosity AGN population (e.g.. 


|Bassani et al. 1 1 999] |Risaliti et al.|1999 |LaMassa et al.|20fT] i. 

Indirect absorption diagnostics predict a larger CT fraction 
for z < 0.5 QS02s than our /V(/.ST4/(-informed distri¬ 
bution. The apparent discrepancy may well be due to in¬ 
direct diagnostics overpredicting the number of CT AGNs. 
Another reconciling factor could be that the quasars unob¬ 
served/undetected with NuSTAR, in general, suffer even heav¬ 
ier absorption than our detected objects. Deeper observations 
at both low (e.g., with Athena', Nandra et al.||2013|l a nd high 
(e.g., with NuSTAR or Astro-H\ Takahash i et al.|2012) X-ray 


energies are needed to reliably distinguish between the above 
scenarios, and thus achieve tighter constraints on /ct for the 
quasar population. 


6 . SUMMARY 


Sensitive high-energy (> 10 keV) NuSTAR observations 
of five optically selected candidate CTQS02s have been pre¬ 
sented, along with broad-band X-ray spectral and multi wave¬ 
length analyses. Similar studies for a further four such objects 
have already been presented in the literature (F14; G14). The 
overall sample of nine z < 0.5 candidate CTQS02s was se¬ 
lected primarily on the basis of multiwavelength evidence for 
absorption by CT (A/j > 1,5 x 10 24 cm -2 ) material along the 
line-of-sight (see Section [A2] i. Our results are summarized as 
follows: 


• Of the five recently observed objects, two are unde¬ 
tected by NuSTAR at 8-24 keV (SDSS J0758+3923 
and 0840+3838), one is weakly detected (net source 
counts S S -24 keV = 38.1/^;?; SDSS J1713+5729), 
and two are strongly detected (Ss -24 keV iZ 90; 
SDSS J1218+4706 and 1243-0232). These represent 
the first detectio ns of these sources at high X-ray ener¬ 
gies; see Section [3~i~| 


• For the two strongly detected targets, spectral modeling 
of the NuSTAR plus archival soft X-ray data suggests 
that the primary transmitted AGN continua are sup¬ 
pressed by extreme levels of photo elec tric absorption. 


with A r H > 10 24 cm , see Section 4.1 For the bright¬ 


est source, SDSS J1218+4706, the relatively high qual¬ 
ity spectral analysis using physically motivated models 
provides strong evidence for CT absorption, wi th a co n- 
tribution from Compton reflection; see Section [4. 1.1 


• For SDSS J1713+5729, the NuSTAR detection likely 
represents the first identification of directl y tran smit- 
ted emission from the AGN; see Section |4.1.3| We 
characterize the 3-24 keV spectrum using the NuS¬ 
TAR band ratio (BRn u ) and estimate a high absorbing 
colu mn density of A/j > 5 x 10 23 cm -2 ; see Sec¬ 
tion |4.2| Notably, the observed 10-40 keV luminos¬ 
ity appears to be extremely suppressed, by a factor of 
« 35, with respect to the intrinsic luminosity, suggest¬ 
ing An > 10 24 cm -2 if purely due to absorption; see 
Section l4.3l 


• For the non detections, column densities of A/j > 
10 24 cm -2 are inferred by comparing the upper limits 
in observed X-ray luminosity (at rest-frame 2-10 keV 
and 10-40 keV) with the intrinsic luminosities expected 
from the mid-IR emission. The majority of NuSTAR- 
observed candidate CTQS02s have X-ray:mid-IR ra¬ 


tios suggesting CT absorption; see Section 4.3 


• For the five objects in the overall NitSTA /(-observed 
candidate CTQS02 sample that are detected at high en¬ 
ergies, the column densities and intrinsic luminosities 
measured from spectral analyses are factors of ss 2.5- 
1600 and « 10-70 higher, respectively, than soft X-ray 
constraints in the literature; see Section [5~i~! 


• Using direct constraints on absorption for 39 QS02s 
studied at X-ray wavelengths, and assuming that the 
Nu ST A /(-detected QS02s are representative of the 
larger sample with evidence for CT absorption, we 
make a prediction for the A/j distribution of optically 
selected QS02s. The distribution is highly skewed to¬ 
ward large column densities (TVh > 10 23 cm -2 ) and 
the predicted CT fraction of /ct = 36^2 % is broadly 
consistent with CXB models. A higher fraction of up 
to 76% is possible if indirect absorp tion diagnostics are 
assumed to be accurate; see Section [5721 


ACKNOWLEDGEMENTS 

We thank the referee for their careful review, which 
has improved this work. We acknowledge financial sup¬ 
port from: the Science and Technology Facilities Coun¬ 
cil (STFC) grants ST/K501979/1 (G.B.F.), ST/J003697/1 
(P.G.), ST/1001573/1 (D.M.A. and A.D.M.); the Feverhulme 
Trust (D.M.A.); Gemini-CONICYT grant 32120009 (R.J.A.); 
the ERC Advanced Grant FEEDBACK at the University of 
Cambridge (J.A.); NSF AST award 1008067 (D.R.B.); the 
NASA Earth and Space Science Fellowship Program, grant 
NNX14AQ07H (M.B.); CONICYT-Chile grants Basal-CATA 
PFB-06/2007 (F.E.B.), FONDECYT 1141218 (F.E.B.), and 
“EMBIGGEN” Anillo ACT1101 (F.E.B.); the Ministry of 
Economy, Development, and Tourism’s Millennium Science 


























17 


Initiative grant IC120009, awarded to The Millennium Insti¬ 
tute of Astrophysics, MAS (F.E.B.); Caltech NuSTAR sub¬ 
contract 44A-1092750 (W.N.B. and B.L.); NASA ADP grant 
NNX10AC99G (W.N.B. and B.L.); the Caltech Kingsley visi¬ 
tor program (A.C.); ASI/INAF grant 1/037/12/0011/13 (A.C., 
S.P., C.V.); NASA ADAP award NNX12AE38G (R.C.H.); 
National Science Foundation grant 1211096 (R.C.H.); and 
Swiss National Science Foundation grant PP00P2_138979/1 
(M.K.). We thank Andrew Ptak and Jianjun Jia for the use¬ 
ful correspondence. This work was supported under NASA 
Contract No. NNG08FD60C, and made use of data from the 


NuSTAR mission, a project led by the California Institute of 
Technology, managed by the Jet Propulsion Laboratory, and 
funded by the National Aeronautics and Space Administra¬ 
tion. We thank the NuSTAR Operations, Software and Cali¬ 
bration teams for support with the execution and analysis of 
these observations. This research has made use of the NuS¬ 
TAR Data Analysis Software (NuSTARDAS) jointly devel¬ 
oped by the ASI Science Data Center (ASDC, Italy) and the 
California Institute of Technology (USA). 

Facilities: Chandra, IRAS, NuSTAR, SDSS, Spitzer, Swift, 
WISE, XMM-Newton. 


REFERENCES 


Akylas, A., Georgakakis, A., Georgantopoulos, I.. Brightman, M., & 

Nandra, K. 2012, A&A, 546, A9S 

Alexander, D. M., Bauer, F. E., Chapman, S. C., et al. 2005, ApJ, 632, 736 
Alexander, D. M., Chary, R.-R.. Pope, A., et al. 2008, ApJ, 687, 835 
Alexander, D. M., Stern, D., Del Moro, A., et al. 2013, ApJ, 773. 125 
Antonucci, R. 1993, ARA&A, 31, 473 

Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference 
Series, Vol. 101, Astronomical Data Analysis Software and Systems V, 
ed. G. H. Jacoby & J. Barnes, 17 

Assef, R. J., Kochanek, C. S., Brodwin, M., et al. 2008. ApJ, 676, 286 
—. 2010, ApJ, 713, 970 

Assef, R. J., Stem, D., Kochanek, C. S., et al. 2013, ApJ, 772, 26 
Assef, R. J., Eisenhardt, P. R. M., Stem, D., et al. 2014, ArXiv e-prints, 
arXiv: 1408.1092 

Balokovic, M., Comastri, A., Harrison, F. A., et al. 2014, ApJ, 794, 111 
Bassani, L., Dadina, M., Maiolino, R., et al. 1999, ApJS, 121, 473 
Bauer, F. E., Yan, L., Sajina, A., & Alexander, D. M. 2010, ApJ, 710, 212 
Bauer, F. E„ Arevalo, P, Walton, D. J., et al. 2014, ArXiv e-prints, 
arXiv: 1411.0670 

Baumgartner, W. H„ Tueller. J., Markwardt, C. B., et al. 2013, ApJS, 207, 19 
Brightman, M„ & Nandra, K. 2011, MNRAS, 413, 1206 
Buchner, J.. Georgakakis, A., Nandra, K., et al. 2015, ApJ, 802, 89 
Burlon, D., Ajello, M., Greiner, J., et al. 2011, ApJ, 728, 58 
Cash, W. 1979. ApJ. 228, 939 

Comastri, A. 2004, in Astrophysics and Space Science Library. Vol. 308, 
Supermassive Black Holes in the Distant Universe, ed. A. J. Barger, 245 
Comastri, A., Setti, G., Zamorani, G., & Hasinger, G. 1995, A&A, 296, 1 
Comastri, A., Ranalli, P, Iwasawa, K., et al. 2011. A&A, 526, L9 
Del Moro, A., Mullaney, J. R., Alexander, D. M., et al. 2014, ApJ, 786, 16 
Donoso, E., Yan. L., Stern, D., & Assef, R. J. 2013, ArXiv e-prints, 
arXiv:1309.2277 

Draper, A. R., & Ballantyne, D. R. 2010, ApJ, 715, L99 
Fabian, A. C. 1999, MNRAS, 308, L39 

Fabian, A. C.. Vasudevan, R. V., Mushotzky, R. F., Winter, L. M„ & 
Reynolds, C. S. 2009, MNRAS, 394, L89 
Fiore, F., Puccetti, S., Brusa, M„ et al. 2009, ApJ, 693, 447 
Gallagher, S. C., Brandt, W. N„ Laor, A., et al. 2001, ApJ, 546, 795 
Gandhi, R, Horst, H., Smette, A., et al. 2009, A&A, 502, 457 
Gandhi, R. Yamada, S„ Ricci. C., et al. 2015, MNRAS, 449, 1845 
Gandhi, P., Lansbury, G. B., Alexander, D. M., et al. 2014, ApJ, 792, 117 
George, I. M.. & Fabian, A. C. 1991, MNRAS, 249, 352 
Ghisellini, G„ Haardt, F„ & Matt, G. 1994, MNRAS, 267, 743 
Gilli, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Gilli, R„ Su, J., Norman, C„ et al. 2011, ApJ, 730, L28 
Goulding. A. D., Alexander. D. M., Bauer, F. E., et al. 2012, ApJ, 755, 5 
Goulding, A. D., Alexander, D. M., Mullaney, J. R., et al. 2011, MNRAS, 
411,1231 

Hainline, K. N., Hickox, R., Greene, J. E., Myers, A. D., & Zakamska, N. L. 
2013, ApJ, 774, 145 

Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103 
Hasinger, G., Miyaji, T., & Schmidt, M. 2005, A&A, 441, 417 
Hazard, C., Mackey, M. B., & Shimmins, A. J. 1963, Nature. 197, 1037 
Heckman, T. M., Ptak, A., Homschemeier, A., & Kauffmann, G. 2005, ApJ, 
634, 161 

Hickox, R. C., Jones, C., Forman, W. R., et al. 2007, ApJ, 671, 1365 
Hopkins, P. F., Hernquist, L., Cox, T. J., & Keres, D. 2008, ApJS, 175, 356 
Humphrey, A., Villar-Martin, M., Sanchez, S. F., et al. 2010, MNRAS, 408, 
LI 

Iwasawa, K., Gilli, R., Vignali, C., et al. 2012, A&A, 546, A84 

Jia, J., Ptak, A.. Heckman. T., & Zakamska, N. L. 2013, ApJ, 777, 27 

Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 


Koss, M., Mushotzky, R., Baumgartner, W., et al. 2013, ApJ, 765, L26 
Kraft, R. P„ Burrows, D. N„ & Nousek, J. A. 1991, ApJ, 374, 344 
Lacy, M.. Ridgway, S. E., Sajina, A., et al. 2015, ArXiv e-prints, 
arXiv:1501.04118 

LaMassa, S. M., Heckman, T. M., & Ptak, A. 2012, ApJ, 758, 82 
LaMassa, S. M., Heckman. T. M., Ptak, A., et al. 2009. ApJ. 705, 568 
—. 2010, ApJ, 720, 786 
—. 2011. ApJ, 729, 52 

LaMassa, S. M„ Yaqoob, T., Ptak, A. F., et al. 2014, ApJ, 787, 61 
Lansbury, G. B., Alexander, D. M., Del Moro, A., et al. 2014, ApJ. 785, 17 
Lanzuisi, G., Piconcelli, E., Fiore, F., et al. 2009, A&A, 498, 67 
Lanzuisi, G., Ranalli, R, Georgantopoulos, I., et al. 2015a, A&A, 573, A137 
Lanzuisi, G., Perna, M., Delvecchio, I., et al. 2015b, ArXiv e-prints, 
arXiv:1505.01153 

Lonsdale Persson, C. J., & Helou, G. 1987, ApJ, 314, 513 

Luo, B„ Brandt, W. N„ Alexander, D. M„ et al. 2013, ApJ, 772, 153 

—. 2014, ApJ, 794, 70 

Lutz, D., Maiolino, R., Spoon, H. W. W., & Moorwood, A. F. M. 2004, 
A&A, 418, 465 

Madsen, K. K., Harrison, F. A., Markwardt, C., et al. 2015, ArXiv e-prints, 
arXiv: 1504.01672 

Magdziarz, R, & Zdziarski, A. A. 1995, MNRAS, 273, 837 
Maiolino, R., Salvati, M.. Bassani, L„ et al. 1998, A&A, 338, 781 
Mateos, S., Alonso-Herrero, A., Carrera, F. J., et al. 2013, MNRAS, 434, 
941 

Mateos, S., Carrera, F. J., Page, M. J., et al. 2010, A&A, 510, A35 
Mateos, S., Carrera, F. J., Alonso-Herrero, A., et al. 2015, ArXiv e-prints, 
arXiv:1501.04335 

Murphy, K. D„ & Yaqoob, T. 2009, MNRAS, 397, 1549 
Nandra, K., Barret, D., Barcons, X., et al. 2013, ArXiv e-prints, 
arXiv: 1306.2307 

Nousek, J. A., & Shue, D. R. 1989, ApJ, 342, 1207 
Ptak, A., Zakamska, N. L., Strauss, M. A., et al. 2006, ApJ, 637, 147 
Reyes, R., Zakamska, N. L., Strauss, M. A., et al. 2008, AJ, 136, 2373 
Risaliti, G. 2002, A&A, 386, 379 

Risaliti, G., Maiolino. R.. & Salvati, M. 1999, ApJ, 522, 157 

Sanders, D. B., Soifer, B. T., Elias, J. H., et al. 1988, ApJ, 325, 74 

Sargsyan, L., Weedman, D., Lebouteiller, V., et al. 2011, ApJ, 730, 19 

Schmidt, M. 1963, Nature, 197, 1040 

Schmidt, M„ & Green, R. F. 1983, ApJ, 269, 352 

Scott, A. E., Stewart, G. C., Mateos, S., et al. 2011, MNRAS, 417, 992 

Simpson, C. 2005, MNRAS, 360, 565 

Smith, R. K., Brickhouse, N. S., Liedahl, D. A., & Raymond, J. C. 2001, 
ApJ, 556, L91 

Stem, D. 2015, ArXiv e-prints, arXiv:1506.00162 
Stem, D., Assef, R. J., Benford, D. J., et al. 2012, ApJ, 753, 30 
Stem, D., Lansbury, G. B., Assef, R. J., et al. 2014. ApJ, 794, 102 
Takahashi, T., Mitsuda, K., Kelley, R., et al. 2012, in Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 
8443, Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, 1 

Teng, S. H„ Brandt, W. N„ Harrison, F. A., et al. 2014. ApJ, 785, 19 
Treister, E„ & Urry, C. M. 2005, ApJ, 630, 115 

Treister, E., Urry, C. M., Schawinski, K., Cardamone, C. N., & Sanders, 

D. B. 2010, ApJ, 722, L238 

Treister, E., Urry, C. M., & Virani, S. 2009, ApJ, 696, 110 
Ueda, Y., Akiyama, M., Hasinger, G., Miyaji, T., & Watson, M. G. 2014, 
ApJ, 786, 104 

Ueda, Y„ Akiyama, M„ Ohta, K„ & Miyaji, T. 2003, ApJ, 598, 886 

Urry, C. M„ & Padovani, P. 1995, PASP, 107, 803 

Vignali, C., Alexander, D. M., & Comastri, A. 2006, MNRAS, 373, 321 


18 


Lansbury et al. 


Vignali, C„ Alexander, D. M„ Gilli, R„ & Pozzi, F. 2010, MNRAS, 404, 48 
Vignali, C„ Mignoli, M„ Gilli, R„ et al. 2014, A&A, 571, A34 
Villar-Martln, M., Cabrera Lavers, A., Bessiere, P., et al. 2012, MNRAS, 
423, 80 

Wachter, K„ Leach, R„ & Kellogg, E. 1979, ApJ, 230, 274 
Weisskopf, M. C., Wu, K„ Trimble, V., et al. 2007, ApJ, 657, 1026 
Werner, M. W„ Roellig, T. L„ Low, F. J., et al. 2004, ApJS, 154, 1 
Wik, D. R., Hornstrup, A., Molendi, S., et al. 2014, ApJ, 792, 48 
Wilkes, B. J., Kuraszkiewicz, J., Haas, M., et al. 2013, ApJ, 773, 15 


Wright, E. L., Eisenhardt, P. R. M.. Mainzer, A. K., et al. 2010, AJ, 140, 
1868 

Wu, J., Brandt, W. N„ Hall, P. B„ et al. 2011, ApJ, 736. 28 
Yaqoob, T. 1997, ApJ, 479, 184 

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579 
Zakamska, N. L., Strauss, M. A., Krolik, J. H„ et al. 2003, AJ, 126, 2125 
—. 2006, AJ, 132, 1496 


APPENDIX 

Here we give further information on the individual NuSTA /^-observed candidate CTQS02s presented in this paper, namely 
relevant multi wavelength properties and features which provide evidence for CT material (Section A.l). In the case of X-ray 
properties, this Section focuses on the low energy (< 10 keV) Chandra and XMM-Newton data which was available prior to 
the NuSTAR observations. For the broad-band X-ray constraints incorporating high energy (> 10 keV) NuSTAR data, which 
generally suggest extreme absorption, see Section Bl In addition, we discuss the identification of a strong Fe Ka line in the 
XMM-Newton spectrum of SDSS J0011+0056, the single Ai/STAR-detection in the exploratory candidate CTQS02 sample pre¬ 
sented by L14 (Section A.2). Lastly, we provide the near-UV to mid-IR photometric data used in the SED modeling (Section A.3). 

A.l. ADDITIONAL INFORMATION FOR INDIVIDUAL OBJECTS 
SDSS J075820.98+392336.0 (z=0.216) 

Fitting an unabsorbed power law model to the XMM-Newton 0.5-10 keV data, we measure a flat effective photon index of 
r eff = 1.1 ± 0.4, indicative of photoelectric absorptio n in the X-ray spectrum. This source is not detected (above the 2.6cr 
confidence level) with NuSTAR at 8-24 keV; see Section |3~T| 


SDSS J084041.08+383819.8 (z=0.313) 

From HST WFPC2 (F814W) imaging (Villar-Mart m et al.|2012 i, the host has a spiral morphology, with evidence for a tidal 
feature. Humphrey et al. (|2010j) included this object in their integral-field observations of six SDSS-QS02s and found spatially 
extended |0 III] and [O II] emission on scales of up to 27 kpc, consistent with being powered by AGN activity (e.g., via shocks or 
radiation). Using the available XMM-Newton 0.5-10 keV data for this object we measure r e ff = 0.7 ±0.1, a low value suggestive 
of heavy absorption. This source is a non detection in the NuSTAR 8-24 keV data; see Section [3TT| 


SDSS J121839.40+470627.7 (z=0.094) 

The 0.5-10 keV XMM-Newton (obsID 0203270201) spectrum is modeled in J13 and LaMassa et al. j( j2012| ), who measure high 
column densities of Ah = 8 .O +4 ® x 10 23 cm -2 and Ah = ( 8 . 7 + 34 ) x 10 23 cm -2 , respectively. J13 measure a strong Fe Ka 
feature at £); ne = 6.4 ± 0.2 keV with EWp e Ka = L7+ 2 'ij keV, consistent with CT absorption. This target is strongly detected 
with NuST AR at 8 -24 keV, allowing relatively detailed, broad-band spectral modeling which extends to high energies (> 10 keV); 
see Section l4.1.1l 


SDSS J124337.34-023200.2 (z=0.281) 


Using HST ACS imaging, Zakamska et al. 
implying an elliptical morphology. 


(j2006j find that the host galaxy light profile is well fit by a de Vaucouleurs profile. 
The host morphology is notably asymmetric. (Zakamska et al. 1(2006) find no evidence for 

AGN" 


and measure a blue excess in the nucleus 


extinction in the host galaxy, suggesting that kpc-scale dust is not obscuring the . 
which may be due to scattering or starburst emission. 

Studying the existing Chandra data, we find an excess of emission at observed-frame sa 5 keV (i.e., rest-frame sa 6.4 keV). 
When fitting the continuum emission with a power law and the excess with a Gaussian component, we measure a rest-frame 


Fe Ka 


= 2.5' 


+4.2 

-2.4 


keV. 


centroid energy compatible with Fe Ka (/i.j iTie = 6.5 _q 2 keV), and a rest-frame equivalent width of EWpe 
Although the emission is consistent with EWp e Ka > 1 keV, which would suggest the presence of CT material, there are too few 
photon counts to rule out low equivalent widths. The object appears to have an extremely flat spectrum, with r e fj = — l.lj^'g for 
the 0.5-8 keV energy band, indicating strong photoelectric absorp tion. This target is strongly detected at 8-24 keV with NuSTAR, 
allowing broad-band X-ray spectral modeling; see Section[4. 1.2 


SDSS J171350.32+572954.9 (z=0.113) 

The mid-IR spectrum, as measured with Spitzer-IRS ( ]Sargsyan et al.|20TT] ), is AGN-dominated and has evide nce for shallow 
silicate (Si) absorption at sa 10 //m. The low energy X-ray properties of this source are detailed in Section |4. 1.3 
an extremely steep spectral shape at 0.5-10 keV (F sa 3) suggests that the weak NuSTAR detection at 
identification of directly transmitted AGN emission from this system. 


To summarise, 
.—24 keV is the first 


A.2. AN IRON LINE IN THE X-RAY SPECTRUM OF SDSS J001111.97+005626.3 

The < 10 keV X-ray spectrum of SDSS J0011+0056 was first presented in J13. L14 extended the X-ray analysis to high 
energies and used the NuSTAR/XMM-Newton band ratio to identify heavy, close to CT, absorption (An ~ 8 x 10 23 cm -2 ). 
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TABLE 5 

Near-Ultraviolet to Mid-Infrared Source Properties 


SDSS I 

0758+3923 

0840+3838 

1034+6001 

1218+4706 

1243-0232 

1713+5729 

u (0.355 /im)“ 

18.967 

± 

0.025 

20.349 

± 

0.179 

16.139 

± 

0.008 

18.727 

± 

0.030 

20.604 

± 

0.116 

18.721 

± 

0.025 

g (0.468 /im)“ 

18.423 

± 

0.008 

19.166 

± 

0.023 

14.743 

± 

0.002 

17.562 

± 

0.008 

19.334 

± 

0.018 

17.480 

± 

0.006 

r (0.616 /im)“ 

17.792 

± 

0.007 

18.021 

± 

0.014 

14.342 

± 

0.002 

16.843 

± 

0.008 

18.015 

± 

0.010 

16.629 

± 

0.004 

i (0.748 (im)“ 

17.629 

± 

0.008 

17.627 

± 

0.013 

13.871 

± 

0.002 

16.386 

± 

0.008 

17.782 

± 

0.012 

16.133 

± 

0.004 

z (0.892 pm) a 

17.706 

± 

0.019 

17.171 

± 

0.026 

13.698 

± 

0.004 

16.180 

± 

0.014 

17.391 

± 

0.029 

16.093 

± 

0.009 

WISE (3.4 fim) b 

13.847 

± 

0.028 

14.322 

± 

0.029 

11.187 

± 

0.024 

12.592 

± 

0.023 

14.762 

± 

0.040 

12.466 

± 

0.023 

WISE (4.6 gm) 6 

12.267 

± 

0.024 

13.549 

± 

0.035 

10.016 

± 

0.021 

11.448 

± 

0.021 

14.348 

± 

0.063 

11.060 

± 

0.021 

WISE (12 pm) b 

8.659: 

± 

0.022 

10.013 

± 

0.041 

6.295: 

± 

0.014 

8.242 : 

± 

0.019 

11.270 

± 

0.155 

7.242 : 


0.015 

Spitzer (3.6 pm) c 













279.10C 

i± 

: 3.333 




Spitzer (4.5 gmf 













258.60C 

' + 

: 3.668 




Spitzer (5.8 /im) c 













280.000 

± 

10.640 




Spitzer (8.0 gmf 













535.300 

± 

14.130 





Notes. 

a SDSS DR7 model magnitudes in the AB sinh system, corrected for Galactic extinction. 

b WISE magnitudes in the Vega system. We use the gmag magnitude for SDSS J1713+5729, and profile-fit magnitudes for the remainder. 
c Spitzer 3.8” diameter aperture flux densities in units of //.ly. 

L14 did not perform detailed spectral modeling, due to the low source counts (« 25 net source counts). However, studying the 
XMM-Newton 0.5-10 keV spectrum we find evidence for an excess at observed-frame « 4.5 keV (i.e., rest-frame ss 6.4 keV). 
Modeling the continuum emission with a power law and the excess with a Gaussian component, the rest-frame line centroid 
energy is in good agreement with that expected for Fe Ka line emission (/7| ilie = 6.4 ±0.1 keV), and the rest-frame equivalent 
width is large (KVVfc Ka = 2.9 ±2% keV). This strong Fe Ka emission suggests CT absorption, and it adds confidence to the high 
column density measured by L14. 


A.3. NEAR ULTRAVIOLET TO MID-INFRARED PHOTOMETRY 

In Table[5]we provide the near-UV to mid-IR photometric data set for the five NuSTA /f-observed QS02s pre sent ed in this work, 
and the one presented in G14 (SDSS J1034+6001). This data set is adopted for the SED modeling in Section [±3] 






